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Abstract—The development of new artificial materials 
is a key issue of the modern world moving toward the 5G 
communication standard and Industry 4.0. Therefore, 
the characterization of these materials (determination of 
their key parameters) plays an important role. The 
dielectric and magnetic constants and their anisotropy 
can be measured in the microwave range; as integral 
characteristics, they give useful information for the whole 
sample and we have shown how they can depend on the 
sample composition, structure, inclusion orientation, 
fabrication technology, conditions for the sample growth 
and forming, etc. We present in this paper a summary of 
our concepts, models, measurement methodologies and 
measurement methods for extraction of the dielectric and 
magnetic parameters of different artificial materials, as a 
result of our experience in the material characterization 
in the Microwave Laboratory of the Faculty of Physics in 
Sofia University, Bulgaria. A set of useful examples is 
presented for the determination of the dielectric and 
magnetic parameters of big groups of artificial materials: 
reinforced substrates, 3D printed dielectrics, ceramics, 
multilayer radomes, foams, absorbers, gradient 
dielectrics and magneto-dielectrics, textile fabrics, meta- 
materials, carbon-containing materials, fresh plant 
tissues, plasmas, etc. A special circumstance in this 
research is the anisotropy of the considered materials 
(different permittivity and/or permeability in different 
directions), which is a very informative parameter and 
gives valuable additional information for the samples. We 
discuss in this paper the origin of the anisotropy of 
variety of materials, how to measure this parameter, how 
to use the anisotropy data in the design and simulations 
of microwave devices and antennas and how the 
anisotropy influences their main characteristics. 


Index Terms—3D-printed dielectrics, absorbers, 
antenna radomes, artificial dielectrics, carbon-containing 
materials, ceramics, conductivity, composites, crystals, 
dielectric anisotropy, ferrites, metamaterials, microwave 
measurement methods, multiferoics, nano-materials, 
permeability, permittivity, substrates, textile fabrics 


I. INTRODUCTION 


The new materials are a key issue of today's world. The 
mechanical and chemical engineers together with the 
physicists can construct now a variety of new artificial 
materials with quite specific properties, never met before in 
nature. As we mentioned in the title, the modern materials 
became reinforced, composite, 3D printed, mixed, nano- 
sized and meta-materials. First of all, these new materials 
will play a very important role for the realization of the 
current trends for automation and data exchange in 
manufacturing technologies within the modular structured 
smart factories or so-called fourth industrial revolution 
Industry 4.0. The second field of the modern materials’ 
implementation is the forthcoming 5G/6G communication 
standards, which will ensure the ability of machines, devices, 
sensors and people to connect and communicate with each 
other via the Internet of Things (IoT) or Internet of People 
(IoP). The third important area of application of the new 
materials is the sustainable development of human societies. 
All these future processes in the development of human 
civilization explain the needs to develop materials with a 
variety of new properties and the needs to have reliable ways 
to characterize these unique properties. 


The characterization of materials usually means to specify 
their properties by different chemical, mechanical, optical, 
spectroscopic, electric, magnetic and other methods. The 
determination of the electromagnetic (EM) properties is also 
important when the new materials can be used in the modern 
electronics, wired and wireless communications, sensor 
networks, for specific EM wave propagation, EM 
compatibility, optical and microwave imaging, etc. The EM 
properties of materials are connected mainly with their 
complex dielectric and magnetic constants and related 
parameters as conductivity, refractive index, wavelength, 
propagation constants, phase/group velocities of waves, 
attenuations, penetration depths, cut-off and threshold 
conditions, etc. in the EM structures, where these materials 
have been incorporated. Due to these reasons, today many 
materials, known with some of their traditional applications 
can be considered as electrodynamic media, e.g. textile 
fabrics as antenna substrates, bio-tissues as electronic 
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components, semiconductor structures as sensors for 
different macroscopic parameters of the environment, etc. 


Our laboratory for Microwave material characterization 
has been orientated to determination mainly of the dielectric 
(and magnetic) properties of the modern artificial materials: 
3D printed dielectrics, reinforced substrates, ceramics, multi- 
layer antenna radomes, foams, absorbers, gradient dielectrics 
and magneto-dielectrics, wearable textile fabrics, metamate- 
rials, carbon-containing materials, liquids, fresh plant tissues, 
plasmas, ferrites, etc. We concentrated our research to deter- 
mination of their dielectric anisotropy (different permittivity 
along different directions), which turned out to be a very 
informative parameter and could be bound with the sample 
composition, structure, inclusions’ orientation, fabrication 
technology, conditions for sample growth and forming, etc. 


In this review article originally published in [1] in two 
parts, we try to summarize the common issues of our many 
years’ research work for characterization of different 
materials in the whole microwave range, including the mm- 
wavelength range. A poster with the same title was presented 
on the interactive University poster exhibition of the Asia- 
Pacific Microwave Conference APMC'2018 [2]. 


The article is conditionally divided into two parts. First of 
all, we present here a classification of the modern materials 
from an EM point of view, their main EM parameters and 
how these parameters could be bounded with the structure, 
inclusions and used fabrication technology. Especially, we 
discuss in the first part the origin of the different types of 
artificial anisotropy with examples. Then we describe the 
developed in the laboratory variety of measurement methods 
for characterization of the material parameters in the 
microwave range (0.5-40 GHz), including several authorship 
methods. We include also the proposed numerical models for 
reliable extraction of the material parameters from the 
obtained measurements results and discuss the basic 
principles and specific techniques for utilization of the 
commercial 3D EM simulators as auxiliary tools for 
extraction of the material parameters of the samples. In the 
second part of this paper, we illustrate our rich experience in 
the area of material characterization by presenting a lot of 
published and new concrete examples for determination of 
the parameters and specific properties of many commercial 
and unique artificial materials: reinforced substrates, 
ceramics and other crystal samples, multi-layer radomes, 
foams, absorbers, gradient dielectrics and magneto- 
dielectrics, textile fabrics, 3D printed dielectrics, metamate- 
rials, carbon-containing materials, nano-composites, fresh 
plant tissues, plasma, ferrites, etc. Finally, we unify the 
common issues of these examples and present conclusions. 


II. MAIN ELECTROMAGNETIC PARAMETERS OF MATERIALS 
IN THE MICROWAVE RANGE AND How THEY CAN BE USED 
FOR MATERIAL CHARACTERIZATION AND TECHNOLOGY 
CONTROL 


A. Classification of Modern Materials 


From the electrodynamic point of view, the materials 
could be characterized on the base of their response on the 
presence of EM fields in each considered electrodynamic 


system, where these materials have been incorporated. This 
response is determined by the properties of the materials, 
described by defining macroscopic scalar complex 
parameters: dielectric constant (permittivity) € = && = (€', — 
jE" )& = &e'(1 — jtano,) and magnetic constant (permea- 
bility) “= to = (ui — Ju") Mo = Molt’C1 — jtano,) of these 
materials (&, Mo — absolute material constants of the vacuum; 
tan6,; tano, — dielectric and magnetic loss tangents; details in 
the next section). The real parts of the relative parameters 
(é'.; u’) can have as positive, as well as negative values for 
different materials. Exactly this circumstance is the basis for 
the nowadays classification of modern materials [3, 4]. First 
of all, most of the standard dielectrics and some magneto- 
dielectrics fall in the group of double-positive (DPS or right- 
hand) materials, for which ¢'’, > 0; wu’. > 0. They still are the 
most occurring electrodynamic media in modern electronics. 
The other two media are also well-known from nature: a 
group of the “epsilon“-negative (ENG) materials — ¢'. < 0; uw‘, 
> 0 and a group of the “mu“-negative (MNG) materials — ¢’, 
> 0; uw’. < 0. Many plasma media (ENG) and ferrites (MNG) 
at certain frequency intervals in external dc magnetic fields 
fall in the last two groups. The fourth group represents the 
double-negative (DNG or left-hand) materials — ¢'. < 0; uw’. < 
0; this is the class of the most rapidly developing engineering 
artificial materials and negative reflective index meta- 
materials with very specific and extra-ordinary properties, 
not existing in nature, which gives us the right to call them 
21“ century electromagnetic and photonic materials, as it is 
pointed in [5]. 


B. Dielectric and Magnetic Properties of Materials 


The complex permittivity €, of the materials in the micro- 
wave region varies between the static (dc) value & at very 
low frequencies (including w = 2nf = 0) up to the optical 
value & at very high frequencies according to the empirical 
Cole-Cole model [6], often used to describe the dielectric 
relaxation in the water and polymers 
a “de Foo (1) 

1+(jor) | 

Exactly due to this reason the “microwave” complex 
permittivity yw = & of the materials could be enough 
informative parameter for the relaxation mechanism 
understanding, specific sample structure, used inclusions and 
other properties of the materials (wis the angular frequency, 
tT — relaxation time; @, = 1/7 -— relaxation frequency; @ is a 
parameter between O and 1; œ = O describes pure Debye 
relaxation model; a> 0 — stretched relaxation model). 


From a more common point of view, the “permittivity ¢ is 
a measure of how an electric field affects, and is affected by 
a dielectric medium, and is determined by the ability of a 
material to electrically polarize in response to the field, and 
thereby reduce the total electric field inside the material” [7]. 
It is directly related to electric susceptibility Ye, which is a 
measure of how easily a dielectric polarizes in response to an 
electric field. They are related to each other through expres- 
sion €, = (1+ y,), where é&, is a scalar quantity in the simple 


case of linear, homogeneous and isotropic materials, while 
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for anisotropic medium, the permittivity is a tensor ¢€.. 


By similar considerations, the permeability u is the degree 
of magnetization that a material obtains in response to an 
applied external magnetic field. The relative permeability 

.=(+y,,) is related to the magnetic susceptibility Xm, 


which is sometimes called volumetric or bulk susceptibility, 
to distinguish it from y, (magnetic mass or specific 
susceptibility) and yy (molar or molar mass susceptibility). 
Again, 4, 1S a scalar quantity for linear, homogeneous and 
isotropic magnetic materials, while when the medium is 
anisotropic (gyrotropic), the permeability is a tensor jZ. . 


In general, as the permittivity (œ), as well as the 
permeability 4(@), are not constants; they both can vary 
with the local position in the medium, sample form 
(depolarization, demagnetization), the frequency of the fields 
applied (electric/magnetic), humidity, temperature, and other 
parameters. In the non-linear media, they may depend on the 
strength of the electric/magnetic field. 


Most of the modern artificial materials could be 
considered as mixtures between two or more isotropic 
components. In this case, very important is the 
homogenization of the resultant material, which can be 
described by effective (we prefer to use the term 
“equivalent”) dielectric or magnetic constant (& ; Heg). The 
process of homogenization is possible to a certain extent 
degree [8] if the characteristic length of the non-resonant 
inclusions in the mixture is sufficiently smaller than the 
wavelength 4 at the operating frequency (typically less than 
A/10 to A/1000). There exist many different technologies to 
mix two dielectrics: series, parallel, series-parallel; layered, 
dispersed (uniform or random), pillar, impregnated, foamed, 
woven, knitted, by additive or subtractive technology, etc. 
The question is how to calculate the resultant effective/ 
equivalent parameters of the mixed material. A set of 
empirical formulas exist [9]; one of the earlier 
approximations is the well-known Maxwell-Garnett 
expression [7], when the mixture is got from small spherical 
dielectric inclusions of permittivity ¢2, embedded in a host 
environment medium of permittivity €. 


E2 — E] 
9 
-V (£2 — £) 


where V> is the normalized volume fractions of the 
inclusions; the total volume is V; + V, = 1. More close to the 
reinforced microwave substrates, textile fabrics, 3D printed 
and layered composites is the Wiener formula for laminar 
mixtures [10] 


(2) 


Eo + 2E] 


1 
_% v, 
E€ +U 











O<u<so, (3) 


Es, TU 


a & +u 


where the parameter u depends on the orientation between 
the layers, namely: 


Le on 
fog E2 ĉl 
Eeg = EnV) + &V); U = 0 (3.2); u = A4 E26] (3.3) 


for series layers (Reuss bound) (3.1); for parallel layers 
(Voigt bound) (3.2) and for isotropic mixtures with randomly 
arranged grains (Bruggman formula) (3.3) see the 
illustration curves in Fig. 2.la,b [A52]. 


The effective (or equivalent) media models are very im- 
portant for the characterization of the different mixed dielec- 
trics. They present the difference between the pure parallel 
and pure perpendicular to the E-field direction dielectric 
inclusions. There exist two main ways to realize multi- 
component isotropic (homogenized) dielectrics with well- 
designed dielectric constant: 1) by random mixing (as for the 
foams) or 2) by constituents with similar dielectric parame- 
ters (as for same microwave substrates and woven/ knitted 
textile fabrics). However, in the case of 3D printing, these 
methods are not technologically applicable; here we should 


use symmetrical unit cells for the infilling medium (see IV). 
3.0 
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Fig. 2.1. Minimal and maximal bounds for the resultant equivalent 
dielectric constant (a) and dielectric loss tangent (b) of two mixed dielectrics 


The next step for generalization of the concept for the EM 
material constants is the introducing of resonance effects in 
the permittivity and permeability behaviour. In the dielectric 
materials, the charges are displaced by an incident electric 
field. The amplitude is given by the restoring forces on the 
electrons and depends strongly on the time-variation of the 
field and this dipole (formed by the displaced electrons and 
the remaining ion) will oscillate. For frequencies around a 
certain frequency @ the induced dipole moment might be 
very large and strongly depend on the frequency. This 
increasing interaction (called “resonance”) of the incident 
EM field with the structure will lead to strong scattering and 
absorption. The behaviour of the oscillating dipole can be 
described on the base of a simple driven harmonic oscillator 
and the solution gives the known relation between the 
polarization P and electric field E expressed by the electric 
susceptibility y., namely P=é)yv,E. The solution to the 


differential equation describing forced resonance is governed 
by a complex Lorentzian function. Taking in mind that the 
electric displacement D is related to the polarization density 
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P by D=éE+P=6)(1+7, )E=€,€9E , the effective (or 
equivalent) resonance dielectric constant by the Lorentz 
model is [7, 8]: 

2 
———; #-+, A 
TO = JOY 
in which qa is the natural resonance frequency, @ is the 
operating frequency, @, is the so-called plasma frequency 
and ¥ is the damping factor in [loss/s] (here also N is the 
number of atoms per unit volume, q is the unit charge, m, 1s 
the effective electron mass). This model can be generalized 
for multiple resonances. 


Eeg (0) = ]+ 


By a similar way, a Lorentz model can be introduced for 
permeability [11]: 
Fo° 


lalo) = 14+ — $=, = 4 (5) 
i Oo- jOYn 


where now @no 1s the magnetic resonance frequency (when 
the magnetic element is the widely used for metamaterials’ 
construction split-ring resonator SRR [12]), %n is the 
magnetic damping factor, and F is a parameter that depends 
on the loop gap in the SRR. 


Both expressions (4, 5) show that the materials can have at 
some frequency intervals either positive or negative values 
for the effective/equivalent dielectric and magnetic 
constants, or even double negative values (DNG) — see 
examples on Fig. 2.2a,b. This approach is widely used in 
modern metamaterial design [13]. Due to these reasons, 
some other parameters, which describe the metamaterials, 
can have anomalous behaviour. For example, the refractive 
index n may obtain negative signs for some metal-containing 
metamaterials [8] (for isotropic and linear media the simple 


expression is n =£, 4, ). 


} El O | Heg(@) n 10 


Ne Ime eq) 
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Fig. 2.2. Illustration of possible frequency behaviour of the real effective 
(equivalent) dielectric &, or magnetic “, constants for (a) artificial 
dielectrics and (b) artificial magnetics with resonance behaviour; (c) the real 
(solid blue) and imaginary (dashed red) parts of the relative permittivity as 
functions of frequency of sample material following the Drude model 


For the special case of metals, the Lorentz model 
transforms into the Drude model. In metals, most electrons 
are free because they are not bound to nuclei (the restoring 
force is negligible and the natural frequency a = 0 in (4)); 
therefore they will oscillate harmonically under the influence 
of an incident EM wave and this can be described by an 
oscillating current density J = cE (for linear media), where o 
is the conductivity. Taking in mind the relations between the 
polarization P, current J and electric field E [7], we can 
introduce a more common definition for the relative 


equivalent permittivity (so-called relative generalized 
permittivity ¢,,(@) ): 
Eeq(®) = Eeqg(@) — JOeg(@)/@, (6) 


where the equivalent complex quantity ¢,,(@) measures the 
polarization response due to bound charges, and o,,(@) is the 
equivalent conductivity, taking in mind the response of the 
free charges. Thus, the generalized permittivity in conductors 
has a large imaginary part (the second term in (6)) and this 
means that the losses are large. But this also means that the 
penetration depth 6 of the EM waves into the sample surface 
is small. The fields attenuate rapidly within the conductor 
and the fields are confined to a thin surface layer; this allows 
the definition of surface current and surface impedance, 
often used by the electronics and electrical engineers. Fig. 
2.2c displays the frequency behaviour of the permittivity in 
the Drude model. As can be seen, at low frequencies the 
imaginary part dominates the response. Indeed, the typical 
conductivity behaviour for metals at low frequencies 
Ei(@) >-jo,.(@)/@ is included in the Drude model in 
the limit @ — O, and the limit for the de conductivity is 
Cue = O Ep / y.. The Drude model is often used to describe 


the optical properties of noble metals, from which most 
optical metamaterials are made. 


All the considered approximations, models and expres- 
sions above clearly show that the material constants are well 
determined by the overall structure of the samples and their 
interaction with the EM waves, by the type of components in 
the resultant complex mixtures, by orientation and position 
of the building blocks and applied technology. The specific 
integral character of the equivalent material constants (&,, 
eq) gives a potential to ensure additional information for the 
material properties to the results, already obtained by 
different local methods (e.g. spectroscopy, SEM images, 
scanning electron microscope, etc.) that are widely used for 
their characterization. In this paper, we concentrate our 
efforts exactly over the reliable characterization of the 
material constants (mainly permittivity) in the microwave 
range as an EM representative of the materials in different 
electrodynamic projects. Moreover, the EM response 
depends also on the material anisotropy. 


C. Anisotropy and Its Origins in the Modern Materials 


The EM anisotropy of a given material can be expressed 
as the existence of different dielectric and magnetic constants 
in different directions. The reasons could be quite different. 
The microwave and optical engineers usually connect the 
anisotropic behaviour with the microwave ceramics and 
optical-glass lenses; this is one of the oldest known types of 
material anisotropy, named crystalline anisotropy for single- 
or poly-crystalline materials (optical glasses, ceramics, 
artificial soft and low-temperature co-fired ceramics LTCC, 
liquid crystals, etc.). In principle, these materials are 
homogeneous, but the anisotropy appears due to the 
existence of different crystallographic axes in the lattices and 
the fact that charges oscillate differently along these 
directions. This anisotropy is usually relatively strong. Our 
early investigations were connected exactly with this type of 
anisotropy of ceramic disks, cylinders, rings and rods by 
different microwave methods [Al, A2, A5-A7]. Then, the 
characterization of crystalline and crystalline-like materials 
continued by the two-resonator method [A34, A53, A56]. 
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Another specific group of plasma/ferrite researchers 
identify the original concept for the material anisotropy with 
the electric or magnetic gyrotropy of gaseous or solid-state 
plasmas (gyro-electric) and ferrites (gyro-magnetic), all in 
external dc magnetic biasing fields. These phenomena 
belong to the so-called induced anisotropy; the dielectric and 
magnetic properties of such natural metamaterials have been 
described with non-symmetric tensors of the permittivity and 
permeability (including non-diagonal components, 
controllable by the external dc biasing magnetic field). The 
ferroelectric materials and films can also be associated with 
this group of electrically gyrotropic materials but in the 
external dc electric biasing field. In our early investigations, 
we performed a lot of attempts for difficult characterization 
of the microwave ferrites with different forms: disk, rod, 
cylinder, prism [Al, A2, A5]; we applied successfully these 
results in our models and methods of design of different 
gyrotropic devices — circulators and isolators (see the review 
papers [A3, A4]). Special types of investigations were 
devoted to ferrite and magneto-dielectric thin films [A15, 
A16, A23] and microwave absorbers [A10, A11, A28, A50]. 
Our contribution to the characterization of gaseous plasma 
media was connected with the development of classical and 
new types of hairpin resonance probes for the determination 
of plasma density [A27, A35, A37]. 


Nowadays, a new class of artificial materials appeared, for 
which the measured anisotropy could be considered more as 
an undesired property. First of all, in this group, we can add 
the commercial engineered reinforced substrates with many 
applications in the modern RF and microwave electronics. 
Characterization of the dielectric anisotropy of such popular 
materials is very important for the reliability of the modern 
design of different planar structures on high-frequency 
substrates, especially in the mm-wavelength range. Our 
research is connected mainly with the characterization of 
these commercial materials ((A8-A10, A12, A13, A17, A18, 
A20-A22, A24, A26, A30-A32, A34, A36, A41, A43, A44)). 
In our practice, we measured a lot of reinforced substrates 
from different manufacturers, which anisotropy varies in a 
wide range: from ~1-2 % up to 25 % (see below for this 
parameter) [A31]. The anisotropy of similar structures has 
been caused by the spatial inhomogeneity in different 
directions between the mixed reinforcing fibres net and the 
applied filling. To this group of anisotropic materials we can 
add the textile fabrics [A48] used for wearable antennas, 
multilayer antenna radomes [A14, A119, A25], some 3D 
printed dielectrics [A46, A52, A53] (see sections IV, VI). 


The controllable anisotropy of the engineered meta- 
materials and “bandgap” materials (i.e. controllable dielectric 
constants in different directions) is another new type of 
artificial anisotropy [8, 14], caused by the chosen shape, 
geometry, size, orientation and arrangement of the unit cells, 
which form the “lattice” of the artificial materials. Extremely 
big could be the anisotropy of metamaterials with metal 
inclusions. Therefore, it 1s very similar by origin with the 
anisotropy, caused from the inhomogeneity of the reinforced 
and composite samples, but now this property is well- 
designed and fully desired, because it ensures unusual 
characteristics of some anisotropic metamaterial devices 


[16]: invisible cloaks, wave concentrators and converters, 
superlenses, special sensors, etc. We already started 
characterization of such metamaterials with big anisotropy 
[A47, A51] (see section VII). 


Let’s consider the terminology, related to the anisotropy. 
In the common case, the anisotropic/gyrotropic materials are 
describing by full tensors; the material tensors are non- 
symmetrical for ferrites and plasma [8] but they will not be 
considered in the paper. There exist three cases of so-called 
dielectric anisotropy of the artificial materials, which 
dielectric constants can be described by diagonal tensors, 
namely: biaxial (7.1) and uni-axial anisotropy (7.2) or pure 
isotropy (7.3) (similar expression could be written for the 
magnetic constants; e.g. for ferrites, multiferoics, etc): 


E O 0 

(é.)=| 0 e 0 (7.1); 
0 O & 
E ay 0 0 

(é.)=| 0 £a 0 (7.2); 

0 O- Ei 

E€ 0 0 

(é.)=|0 « o0 (7.3) 
0 0 6, 

whore eo Pea ep og a G Ea 


E, Z En X Ey F Ez 


Typical biaxial anisotropic materials are the crystals, 
optical glasses, liquid crystals, some nanocomposites, many 
metamaterials with non-symmetrical unit cells, 3D printed 
dielectrics, tick textile fabrics, plant and other bio-tissues, 
which can be described with three different scalar dielectric 
constants along the axes Ox, Oy and Oz. The uni-axial aniso- 
tropy has been expressed in relatively thin structures with 
plane symmetry, e.g. materials like the reinforced substrates, 
hard or soft artificial ceramics, LTCC’s, multi-layer antenna 
radomes, some polymers, gradient dielectrics and absorbers. 
The near-to-isotropic behaviour is a property of the well- 
homogenized and foam-like materials and metamaterials: 
foams, bulk plastics, injection-moulded dielectrics, non- 
woven substrates, foamed absorbers, etc. In the next sections 
of this paper, we will consider many concrete examples. 


Finally, the so-called magneto-electric materials appear 
among modern artificial materials [7, 16]. The common issue 
for these materials is the magneto-electric coupling; now, the 


complex electric and magnetic flux densities D and B in 
the so-called bi-anisotropic magneto-electric materials can be 
expressed as 


B=()£+(2)B (8.1) and B=(H)A+(C)E 62 


A 


where (é) : (i) are the permittivity and permeability tensors, 


while (2), è ) are tensors of the corresponding magneto- 


electric coupling coefficients. When the parameters €, u, é ¢ 
are scalars, the materials become bi-isotropic. The above 
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relations show that the equivalent permittivity in such 
materials can depend not only on the response to the E field 
but also on the response to the H field in the considered 
media (v.v. for the equivalent permeability). We already 
managed to characterize in our practice such samples [A51, 
A58-A61], but this research is in the beginning. 


D. Numerical Methods for Estimation of the Dielectric 
Anisotropy of Artificial Materials 


The discussions in the previous section express the 
significance of anisotropy for each complex material, not 
only for metamaterials. Therefore, it is important to may 
predict the degree of anisotropy in different materials, but 
this is not possible in all cases. Full-wave analysis of 
anisotropic structures is possible in a few simple cases. The 
crystalline type of anisotropy can be bound to the lattice 
symmetry — isotropy for crystals with cubic symmetry; uni- 
axial anisotropy — with hexagonal, tetragonal symmetry; bi- 
axial anisotropy — with monoclinic, triclinic, orthorhombic 
symmetry [8]. Another very popular model for crystalline 
anisotropy in high-g ceramics is based on the lumped- 
element approach: the building blocks as electric dipoles can 
be considered at the macroscopic level as a set of parallel or 
series capacitors (for example see the expressions 3.1-2), but 
the calculated anisotropy is usually bigger than the measured 
one. Unfortunately, the other effective-media expressions for 
the resultant permittivity in different mixtures cannot give 
the actual anisotropy [9, 10]. There exist also attempts to 
bind the degree of symmetry of the unit cells for evaluation 
of the all-dielectric 3D printed metamaterials [8, 14]. Of 
course, the numerical simulations by the nowadays 3D EM 
simulators are also wide-spread; many microwave and 
optical engineers successfully apply this approach, and 
therefore, the big software developers pay increasing 
attention the new versions of the commercial 3D simulators 
to may fluently “cope” with the metamaterials. 


In [A46, A48, A52], we proposed a numerical method for 
reliable prediction of the dielectric constants of artificially- 
constructed materials (including all-dielectric metamaterials 
or even samples with metal inclusions) along all three axes in 
the Cartesian system by resonance method in desired fre- 
quency intervals. The idea of this method is to build a well- 
designed unit cell, to reproduce it in a hosting isotropic sub- 
strate and to put the whole sample in a rectangular resonator 
(Fig. 2.3), which supports TE and TM modes with three 
mutually perpendicular directions of the E fields — see illus- 
trations in Fig. 2.4. Independent extraction of the resultant 
dielectric constant along all three axes is possible after 
replacing of the bi- or uni-anisotropic sample under interest 
with an equivalent isotropic sample. We successfully applied 
the described method for many artificial samples: textile fab- 
rics [A47], 3D printed dielectrics [A52, A54], metamaterials 
with surface metal inclusions [A46, A51]. Different unit cells 
could be used for these numerical investigations: spheres, 
cubes, cylinders, prisms, disks and some combinations 
between them, made by isotropic dielectrics or by metals. 
We use for a quantitative measure of the resulting dielectric 
constant/loss tangent anisotropy (parameters AA,, AAtans;) for 
bi-/uni-axial anisotropy the following expressions 





(a) (b) 


Fig. 2.3. (a) Universal unit cell with cylinders; (b) constructed artificial 
sample with repeated unit cells (inclusions) in a hosting isotropic substrate; 
(c) artificial sample in a rectangular resonator 


Oz = E'perp 


E 
E'ox 5 E'oy = E par 
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(a) (b) (c) (d) 
Fig. 2.4. Rectangular resonator with an artificial sample (knitted fabrics), 
which supports different modes with mutually perpendicular E fields along 
the axes: (a) Ox; (b) Oy; (c) Oz and (d) in plane Oxy 
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Fig. 2.5. Dielectric constant anisotropy of artificial nets from simple 
inclusions as regular cylinders, orientated along axes Ox, Oy, Oz or in plane 
Oxy; the E-fields' orientation is fixed either along x (Epar) or along 02 (Eperp) 


AA, ox/y =2(E'ox/ y E02 WE oxiytE 0z ), (9.1) 


MA an 0€,0x/ y = 2(tan 0654 y —tan O- Oz ) /(tan On Ox/ y +tan On Vz ) 
(9.2) 


The dependencies for A,, % versus the ratio &n¢//&s5 In Fig. 
2.) well illustrate a simple but informative case for 
anisotropy prediction: inclusions as dielectric cylinders with 
permittivity &,.; orientated along the axes Ox, Oy or Oz, 
keeping the E fields parallel or perpendicular to the sample 
surface by exiting fixed mode in the rectangular resonator 
(Es. 18 the permittivity of the hosting isotropic medium). The 
simulations show that the resultant permittivity &,,,0, of the 
artificial mixture increases along the cylinders when the E 
field is orientated along the same axis, in comparison to the 
dielectric constant & aroy, perpendicularly to the cylinder 
axis. This effect increases the anisotropy A, 0, of the whole 
artificial substrate along the cylinder axis with the rise of 
ratiO Enci/ Eso, While the anisotropy As oy is close to zero. The 
effect of anisotropy for &nci)} &so Increases for a “fishnet” of 
cylinders. The positive AA, sign means that 4, > E&erp, While 
the negative sign — that S4, < Eerp. The last effect happens, 
when the cylinders are orientated perpendicularly to the 
sample. The presented above investigations have been not 
published jet, however, we confirmed numerically and 
experimentally the described rules for anisotropy of textile 
fabrics [A48], reinforced substrates [A31], multilayer 
antenna radomes [A14, A25], 3D printed dielectrics [A54], 
metamaterials [A46]. Similar investigations are helpful for 
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the design of homogenized and near-to-isotropic 3D printed 
dielectrics for antenna applications [A52]. 


Il]. METHODS FOR EXPERIMENTAL CHARACTERIZATION OF 
DIELECTRIC PARAMETERS OF MATERIALS IN THE 
MICROWAVE RANGE (INCLUDING THEIR ANISOTROPY) 


Nowadays, a lot of measurement methods exists for 
characterization of artificial materials, some of them 
proposed as reference ones (approved from the leading world 
standardization organizations as IPC, NIST, NPL [17-19]), 
or developed by different Universities and scientific research 
organizations. This is a world wealth, nevertheless that some 
methods can give more or less different results for the same 
materials. Three types of methods can be classified in the 
microwave range: free-space, waveguide (transmission-line) 
and resonance methods — Fig. 3.la-d. The first two are 
considered as broadband, but less accurate methods, while 
the resonance methods are specified as more accurate, but 
narrowband ones. A good survey of the microwave methods 
for material characterization has been done in [20]. We 
consider selectively in this paper the promising methods with 
an accent to these, which can measure anisotropy and have 
been already realized in the Microwave Lab in the Faculty of 
Physics of Sofia University “St. Kliment Ohridski” [1, 2]. 


A. Resonance Methods 


Definitely, the resonance methods are considered more 
accurate than the transmission-line methods [21]. The main 
reason is the fact that the sample under interest is a part of a 
bigger resonance area (Fig. 3.1d) and its dielectric and/or 
magnetic constants influence relatively strongly the resonan- 
ce parameters of the excited modes — resonance frequency fo 
and unloaded quality factor Qo of the empty resonator, 
shifting them to f, and Q,, when the sample is placed inside 
(in the case of volume resonators; in other cases the sample 
itself could be a high-Q resonator). The big accuracy is 
achieved since the frequency or frequency differences can be 
measured with smaller errors in the microwave range, 
typically less than 0.01% in the X band. However, other 
factors can decrease the measurement accuracy, e. g. actual 
sample shape and sizes, sample positioning in the resonator, 
environmental factors, etc. Our short survey of the resonance 
methods is concentrated over their ability to measure the 
sample anisotropy, considering the expressions (9.1-9.2). 


A measuring resonator can react to the sample biaxial 
anisotropy if it supports at least three independent modes 
with mutually perpendicular E (or H) fields along the 
Cartesian axes Ox, Oy and Oz in a given frequency interval. A 
good illustration is the resonance perturbation method [22], 
suitable for characterization of the biaxial anisotropy of 
liquid crystal polymers (LCP), which have substantial 
potential for application in the modern microwave 
electronics. The method uses three waveguide modes TEoj, 
TMi? and TEj;9; with mutually perpendicular E-field 
orientations in a non-standard rectangular cavity to extract 
the permittivity values &,, 6, and &,. We apply rectangular 
resonators made by standard waveguides that support 
dominant modes TEjo, (p = 1-6), to extract either dielectric 
or magnetic parameters of microwave absorbers (Fig. 3.2). 


a 


Fig. 3.1. Three popular methods used for characterization of flat samples 
under test (SUT): (a) free-space method in transmission regime; (b) in 
reflection regime; (c) waveguide method; (d) resonance method 





Due to the coarse character of the known perturbation 
expressions [20], we apply more accurate 3D simulations. 


The volume resonators also have good potential to 
determine the sample anisotropy, but there is no universal 
solution. The parallel dielectric parameters &,4/tan ôzpar can 
be measured by TE-mode resonators (as R1 in Fig. 3.3a) 
(classical Courtney’s method [23], Kent’s evanescent-mode 
tester [24], split-cylinder resonator SCR [25], split-post 
dielectric resonator [26], etc.). The perpendicular parameters 
Eperp/taNOzperp Can be estimated by TM-mode resonators (as 
R2) (e.g. as in [A21, 27], reentrant cavities [28]). In fact, 
only a few publications have been directly dedicated to 
anisotropy measurements. Whispering-gallery modes with 
high-azimuthal index in a single dielectric resonator with 
fields far from the screens could be used for anisotropy 
measurement of ultra and extremely low-loss materials [29]. 
The method presented in [30] is based on resonance 
measurements of one prism sample in three mutually 
perpendicular orientations in a rectangular resonator (as in 
[A53, A56] for cylindrical resonators). Very important is the 
problem for reference determination of the proven dielectric 
anisotropy of the commercial microwave substrates, because 
the reference IPC TM-650 2.5.5.5 [17] (clamped stripline 
resonator test method) cannot give this parameter. The 
modified IPC TM-650 Bereskin’s method [31] is difficult for 
realization (many substrates have to be stacked horizontally 
and then a single sample has to be produced by a vertical 
cut); only a few substrate producers apply this method. 
Nowadays, a broadband method is proposed [32], based on 
numerical extraction of the parallel and perpendicular 
dielectric constants from data for the even and odd modes in 
coupled microstrip resonators printed on the substrate under 
test. However, the authors make a wrong conclusion about 
the ability to determine the parallel permittivity of a substrate 
(in fact, they determine the so-called design or microstrip- 
like value and this is proven by our measurements). Similar 
to the Bereskin's method (in stripline geometry) is the 
method, proposed in [33] and implemented by substrate- 
integrated waveguide (SIW) resonators. The novelty here is 
the possibility to increase the working frequency up to mm- 
wavelength range. A common disadvantage of the methods, 
based on the printing of measurement resonators on the 
substrate under test, is the difficult extraction of the dielectric 
loss tangent of raw substrates by taking into account the 
influence of the quality of metallization (the problem is 
discussed in IIIC). Recently, a new TM-mode method 
(balanced-type circular disk resonator [34, 35]) shows big 
potential for characterization of the perpendicular dielectric 
parameters of flat samples up to 110 GHz (now reliable 
results have been obtained in the interval 20-70 GHz). We 
just started to develop this method to help substrate anisotro- 
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Fig. 3.2. Perturbation resonance method for determination of dielectric 
and magnetic parameters of microwave absorbers under test (AUT). Inset: 
simulated E-field distribution of the TEj9s; mode in % part of the resonator 
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Fig. 3.3. (a) TE- (R1; SCR) and TM-mode (R2) cylinder resonators; (b) 
pairs with equal diameters ¢30 mm (R1&R2') or with close resonance frequ- 
encies (R1 &R2); (c) mode spectrum in R1 and R2 resonators (0-25 GHz) 
py determination in the millimetre bands for the 5G standard. 
Combined with the free-space or open-resonator method for 
determination of the parallel dielectric parameters, they both 
can be a good base for characterization of the material 
anisotropy up to mm- and sub-mm wavelength ranges. 


B. Characterization of the Dielectric Anisotropy by the 
Two-Resonator Method and Its Modifications 


In 2000-2002 we first detected the dielectric anisotropy of 
commercial microwave substrates and developed an author- 
ship method for its determination — the two-resonator method 
[A21]. Our first publications [A6, A10] have been based on a 
perturbation approach, but then we developed more accurate 
full-wave analytical models of the measurement resonators 
and for the extraction procedure of the complex dielectric 
constants of uni-axial anisotropic single and multilayer 
samples [A17, A21, A26]. Applying this method we 
managed to accumulate in a short period (2004-2008) a 
tremendous amount of data for the real anisotropy on many 
commercial substrates from different producers (used in our 
antenna projects [A29, A33, A38, A39, A43, A57, A61]. 
However, we have encountered many difficulties for a long 
time in trying to impose the idea of having such a feature 
(anisotropy) of these most spread electronic components 
(microwave substrates) as among the big substrate producers, 
as well as among the most of the users (we discuss this 
“story” for the substrate anisotropy in IVA). 


Originally, the two-resonator method [A21, A31] is based 
on two consecutive measurements with/without samples of 
the resonance characteristics of two symmetrical modes: 
TE; 1n R1 resonator (with E fields parallel to the substrate) 
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and TMoio in R2 resonator (E fields perpendicular to the 
substrate) — Fig. 3.3a and Fig. 3.4a. To achieve the best 
sensitivity of this method the sample is placed in the middle 
of R1, while it lies on the resonator bottom in R2; the chosen 
lower height H> of R2 well separates the TMo;9 and TE), 
resonances. The method can be realized in two options: 1) 
for samples with equal diameters D; = D, (e.g. 30 mm), but 
the TMo;9 mode has lower resonance frequency (7.65 GHz) 
than for TE,o;; mode (13.15 GHz); or 2) for achieving of 
close resonance frequencies (e.g. in the interval 12.6-13.15 
GHz), but applying resonators with different diameters (D; = 
30 mm; D» = 18.1 mm) — see the examples in Fig. 3.3b,c. 
The considered here frequency band of the described method 
can be expanded to higher frequencies by using several 
selected high-order modes: e.g. TEo13, TEo21, TEo31 in R1 and 
TMo20, TEo30 1n R2; however, they should be well-identified 
in the mode spectrum (see the modes marked with arrows in 
Fig. 3.3c). The other way for frequency range widening is to 
increase or decrease the resonator (and sample) diameter. 
The diameter decreasing can increase the operating 
frequencies (e.g. for TEp;; in RI: fı ~ 21.8 GHz at D; = 18.1 
mm; fı ~ 28.7 GHz at D; = 15.0 mm; fi ~ 39.3 GHz at D; = 
10.0 mm; for TMo10 in R2: fa ~ 22.8 GHz at D, = 10.0 mm; fh 
~ 28.5 GHz at D2 = 8.0 mm; by higher-order modes in RI: 
up to 70-80 GHz; in R2: up to 65-70 GHz). However, the 
decreasing of resonator diameter leads to a decreasing of 
measurement accuracy, especially for the dielectric loss 
tangent, due to the smaller volume and lower Q factors in 
both resonators. Conversely, the diameter increasing leads to 
resonance frequency decrease and we managed to determine 
the substrate anisotropy by measurements even at 1.4-2.5 
GHz [A48], but the sample sizes become too big and the 
accuracy again decreases. There is another way to decrease 
the frequency interval for anisotropy determination by small 
samples — to use another pair of R1I/R2 resonators with 
naturally lower resonance frequencies based on split-coaxial 
resonator as R1 [A24] and reentrant resonator as R2 [A22] — 
see Fig. 3.4b. With this new pair of resonators, we again 
decreased the frequencies for determination of substrate 
anisotropy up to 1.7 GHz. In general, we successfully 
widened the applicability of two-resonator method by 
introducing in [A30, A31] the approach to use different pairs 
of measurement resonators from type RI/R2 and to obtain 
unique results for many commercial substrates, which was 
notified in [37, 38]. The third pair with split-dielectric post 
resonators SPDR(e/m) with TE/TM modes with added sub- 


R1: C R2: C R1: SCoaxR R2: ReR 
A 
M n T e 
l a) (b) 
1 D2 
R1: SPDR(e) R2:SPDR(m) EE sample 


metal walls 


ea Ca M dielectric resonator 


(c) 

Fig. 3.4. Two-resonator method, realized by three different pairs of 
TE/TM resonators: (a) cylinder resonators (R1:C or R1:SCR and R2:C); (b) 
split coaxial resonator (R1:SCoaxR) and re-entrant resonator (R2:ReR); (c) 
split post dielectric resonators (R1:SPDR(e) and R2:SPDR(m)) 
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Fig. 3.5. Dependencies of the normalized resonance frequency (a) and 
normalized Q-factor (b) of the dominant mode in the cylinder measurement 
resonators CR1 and CR2 
strate sample allows us preliminary to choose a frequency for 
measurements by the used high-Q dielectric resonators (Fig. 
3.4c). We successfully applied the authorship two-resonator 
method for many different materials in addition to the micro- 
wave substrates, like multilayer composites, thin films, 3D 
printed dielectrics, textile fabrics, and even meta-materials 
(with less success) — details will be discussed later on. 


The measurement errors of the two-resonator method have 
been estimated in [A21] as ~2.5-3% for AA, and ~10-12% 
for AAjtans; in the case of multilayer materials in the Ku band 
(or 1.5% and 8-9% for single-layer substrates). These values 
fix the margins for determination of the practical isotropy of 
materials by the two-resonator method. A big contribution to 
the improvement of measurement accuracy is the concept for 
the equivalent parameters of measuring resonators [A30, 
A31]. Usually, the calculated and measured resonance 
parameters of the empty resonators do not fully coincide, 
foca fomeass Qocatc# Qomeas. The reasons are different: size 
uncertainty, the influence of the coupling loops, tuning 
screws, eccentricity, surface cleanness and roughness, daily 
temperature variations, etc. The problem can be overcome 
namely by the introduction of equivalent resonator 
dimensions and equivalent walls’ conductivity. The 
equivalent geometrical parameters (D or H) could be chosen 
on the base of a simple rule: the variation of which parameter 
influences more the resonant frequency of the empty cavity 
without any additional components? In the case of circular 
resonators the diameter changes “provoke” bigger variations 
of the resonance frequency and Q factor (see dependencies in 
Fig. 3.5). If we tune the selected equivalent parameters in the 
constructed stylized 3D models of empty resonators, we can 
achieve a practical coincidence between the calculated and 
the measured resonance parameters: focale ~ fomeass Qocale ~ 
Qomeas. Thus, in the case of pair of cylinder resonators R1I/R2 
their 3D models include equivalent diameter D,,12 (instead 
Di2), actual height H,,. and equivalent wall conductivity 
Oegi2- Additional equivalent parameters could be introduced 
for the other two pairs of used resonators [A31]. 


C. Free-Space Method and Waveguide (Transmission-Line ) 
Methods 


The measurement setup of the free-space method (in 
transmission or in reflection regimes), consists of two or one 
horn antennas and the flat sample is placed between them or 
between the antenna and a reflecting screen — see Fig. 3.1a,b. 
The measurements of the phase delay and the losses with and 
without sample allow extracting the dielectric constant and 
dielectric loss tangent. This method is very useful for 


measurements in the mm- and sub-mm-wave ranges but 
needs relatively big samples to avoid the diffraction effects 
near to the sample edge. The main problems are the reducing 
of the measurement errors due to multiple reflections and 
sample misalignment. By suitable extraction procedure, these 
errors could be reduced or measurements in a TDR gate 
could be used. Due to the TEM mode propagating in the 
structure, it can be used for extraction only of the parallel 
dielectric parameters of the sample in a transmission regime. 
The method is not well applicable for low-loss materials, but 
it is very popular for characterization of new unknown 
materials, including metamaterials [A51], A58, A59], 
composite materials as antenna radomes [A19, A25, A32] 
and absorbers [A28, A50, A61], all realized in our laboratory 
from X- to Ka-bands. An alternative of the free-space 
method is the resonance quasi-optical (Fabry-Perot) inter- 
ferometer, realized in our laboratory in Ka-band; it can again 
measure the parallel dielectric parameters in mm- and sub- 
mm wavelength ranges even for low-loss materials. 


The waveguide methods are also very popular in the dm-, 
cm-, mm-wavelength ranges and beyond, applying rectangu- 
lar, circular or coaxial waveguides in transmission or reflec- 
tion regimes. For rectangular waveguide in transmission 
regime, the thin sample can be placed between two 
waveguide flanges (as in Fig. 3.1c) or into the waveguide 
fitting its cross-section; the measurement of phase delay and 
insertion losses allows extraction of the sample dielectric 
parameters. Inserted in the waveguide aperture bulk samples 
can act as resonators in the corresponding frequency range; 
the extraction of the dielectric parameters is possible by 
measurement of the resonance frequency and Q factor. For 
thin samples in reflection regime, both dielectric and 
magnetic parameters can be extracted if the sample is placed 
in two fixed positions at a quarter-wavelength distance. 
When the sample is placed in the natural vertical position 
according to the wider wall of the waveguide, the method 
gives the parallel dielectric parameters; if the sample lies 
horizontally on this wall, the method extracts the perpen- 
dicular parameters. Thus, the waveguide method is one of 
the few relatively broadband methods, which can determine 
the anisotropy of the sample. Due to this reason, it is very 
popular for characterization of 3D printed metamaterials 
(e.g. see [14]), including in our laboratory [A28, A51]. 


The waveguide method can be also realized by planar 
transmission lines, printed on the substrate under interest. In 
this variant, it is known as the popular differential phase- 
length (DPL) method and it is one of the most spread 
methods for reliable microwave substrate characterization on 
the base of microstrip lines (MSL). The realization is simple: 
So; parameters of two or more 50-Ohms MSL (or other 
planar lines) with different lengths have to be measured and 
the phase difference AangS>,, deg and difference between the 
insertion losses AS>;, dB have to be determined for each pair 
of lines with length difference AL (Fig. 3.6a). Then the 
important parameters, effective dielectric constant &, and 
attenuation a, can be calculated by the simple expressions: 


2 

AangS,;.c 
ep =| 1a 10.1 
A aa) es 
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a,dB/cm=AS,,,dB/AL,cm. (10.2) 


Fig. 3.6b,c presents an example for measured amplitude 
and phase S,; dependencies for a selected from the market 
commercial substrate (conditionally marked as “Substrate 
3.0”) of thickness 20 mils with catalogue dielectric constant 
3.0. Using totally 6 measured combinations for a set of 4 
MSLs on this substrate with different thicknesses (5, 10, 20 
and 30 mils) useful frequency dependencies for the effective 
dielectric constant and attenuation of MSL on this substrate 
can be obtained by (10.1-10.2) — shown in Fig. 3.7. In 
addition to the DPL method, other planar resonance methods 
are used to draw similar curves [A42]. The obtained 
dependencies are important for the extraction of the 
equivalent dielectric parameters — equivalent dielectric 
constant &, (in this case: MSL value & 4 msz, shown in Fig. 
3.8a) and the corresponding equivalent dielectric loss tangent 
tan 0, g. However, in the last case, the imperfect metallization 
influences the MSL attenuation in a big degree [37]; there- 
fore, the conductor and radiation losses have to be taken into 
account for more accurate extraction of tand;.g [A41, A55]. 


For more than 10 years now, the equivalent dielectric con- 
stant &, is an important parameter for reliable EM 
description of each reinforced substrate and each planar line 
printed on it. This approach allows an useful replacing of this 
anisotropic substrate with an isotropic equivalent practically 
in each RF design project, applying equivalent dielectric 
constant £4 and equivalent dielectric loss tangent tan 0,4 (see 
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Fig. 3.7. Measured dependencies of & us, (a) and attenuation (b) by 
DPL and resonance methods for Substrate 3.0 with different thicknesses 
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Fig. 3.9. Measured & (a) and extracted &, (b) dependencies for six 50- 
Ohms planar transmission lines for Substrate 3.0 (20 mils); (c) averaged by 
the thickness and frequency &, values of six planar lines for Substrate 3.0 
(Dk = 3) and other Substrate 3.38 (Dk = 3.38). Legend: SL — stripline; PS — 
paired strips; MSL — microstrip lines; CPWG — coplanar waveguide; 
GCPWG -— grounded CPWG; SIL — slot line. 


our concept described in [A31, A42]). This is an efficient 
approach and the leading PCB producers and developers of 
the 3D EM simulators started to share with the substrate 
users the needed information for the equivalent dielectric 
constant, often called “design Dk” value and usually 
obtained by the DPL method for MSL. Fig. 3.85 presents a 
useful picture for the considered Substrate 3.0 — a compare- 
son between the measured &4,, Eerp aNd &g musz for different 
thickness. The results clearly show an important fact — the 
producers have managed to decrease the anisotropy for thin 
substrates, which are applicable exactly in the mm-wave- 
length range, where the anisotropy has a bigger influence. 

The main problem is that the parameter ¢, of each 
substrate depends on the type of the printed planar line on 
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this substrate (see dependencies in Fig. 3.9a,b), 1.e. from the 
concrete distribution of the parallel and perpendicular 
electric fields of the dominant propagation mode in a given 
planar line. We managed to prove for the first time this fact 
in our work [A31] and then developed a special concept in 
[A42]. It turned out that this fact is a serious inconvenience 
for the application of the equivalent dielectric constant 
approach (design Dk). However, we presented a partial solu- 
tion based on the introduction of three groups of equivalent 
values — see the information, presented in Fig. 3.9c. The 
proposal is RF designers to use in their planar device projects 
three different values of the equivalent dielectric constant £q 
depending on the dominant planar structure: I — “stripline 
value” Eq SL ~ Eperp (OF process, specification Dk, applicable 
for stripline and SIW); If — “microstrip value” &¢ mst ~ (erp 
+ Eperp)/2 (or design Dk, suitable for microstrip line, coupled 
MSL, paired strips, grounded coplanar waveguide and 
similar microstrip-like lines) and IIIf — “coplanar value” 
Eeq_CPWG S ar (Or a new coplanar design Dk, applicable for 
slot-based planar lines: coplanar waveguide, slot and fin- 
line). The first two values are already in use namely as 
“process” and “design” values [38, 39] (our concept only 
explain the origin of these values), but the third value is a 
new parameter — it is suitable for coplanar and slot-line based 
planar structures, used mostly in the integrated circuits. 
Many new publications confirm the usefulness of this our 
concept for reliable design of planar devices with equivalent 
dielectric constants (more details are discussed in IV). 


D. Utilization of 3D Electromagnetic Simulators as an 
Auxiliary Tool for Material Characterization by 
Resonance Measurements 


The modern material characterization needs the utilization 
of powerful numerical tools for obtaining accurate results 
after modelling of very sophisticated measuring structures. 
Such software instruments could be the 3D EM simulators, 
which demonstrate significant capabilities in the modern RF 
design. Considering recent publications in the area of 
material characterization, it 1s easy to establish that the 3D 
simulators have been successfully applied for measurement 
purposes, too. The ability to use commercial frequency- 
domain 3D simulators as auxiliary tools for accurate 
measurements of the substrate anisotropy based on the two- 
resonator method has been demonstrated for the first time in 
[A19]. Then, this option was developed for all types of 
considered resonators [A31], following few principles — 
simplicity, accuracy and fast simulations. Illustrative 3D 
models for the simplest structures, used in the two-resonator 
method (cylinder resonators R1 and R2), are drawn in Fig. 
3.10a. Three main rules have been accepted to build these 
models for accurate and time-effective processing of the 
measured resonance parameters — 1) a stylized drawing of 
the resonators’ body with equivalent diameters (D,,1,2), actual 
height H,. and equivalent wall conductivity o,¢12; 2) an 
optimized number of line segments for construction of the 
curved cylindrical surfaces (line size less than 1/16 wave- 
length) and 3) a suitable for the operating mode splitting (1/4 
or 1/8 from the whole resonator body; Fig. 3.105), 
accompanied by appropriate boundary conditions at the cut- 
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Fig. 3.10. (a) Equivalent 3D models of two cylinder resonators R1/R2 
and boundary conditions (legend: 1 — finite conductivity (o.,); 2 — E-field 


symmetry; 3 — H-field symmetry); (b) Simulated E-field scalar and vector 
distribution in the measurement resonators (1/8 R1 and 1/4 R2) 
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Fig. 3.11. Dependencies of the normalized resonance frequency (a) and 
Q-factors (b) of the used resonance modes for anisotropic and isotropic 
samples versus dielectric anisotropy AA, , AAtanse 
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off planes. Although the real resonators have the necessary 
coupling elements and holes, the resonator bodies could be 
introduced into the model as pure closed cylinders and this 
approach allows applying the Eigenmode solver of the 
modern 3D simulators. The utilization of the Eigenmode 
option for obtaining of the resonance frequency and the 
unloaded Q factor (notwithstanding that the modelled reso- 
nator is not fully realistic) considerably facilitates the aniso- 
tropy measurement procedure assisted by 3D simulators, if 
equivalent parameters have been additionally introduced and 
symmetrical resonator splitting has been done. 


These optimized models increase the accuracy of the two- 
resonator methods for anisotropy measurements and allow us 
to evaluate the sensitivity of the corresponding resonators 
(with the excited TE or TM modes) to the actual sample 
anisotropy — see dependencies in Fig. 3.11. The analysis of 
the obtained results shows that the chosen modes in the 
different pairs of measurement resonators have very good 
selectivity to the corresponding dielectric constant of the 
sample in the resonator along a given direction with an 
uncertainty less than +0.25 % for the normalized resonance 
frequency faniso/fiso and less than +0.5 % for the normalized Q 
factor Qaniso/Qiso, When the anisotropy varies in the intervals 
MA, AAtansg = £10 % even for relatively thick samples 
(excepting for SCR). 


E. Problems with Reliable 3D Simulations of Microstrip 
Lines on Anisotropic Substrates 


The microwave engineers have two possibilities to take 
into account the material anisotropy in their numerical design 
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projects — to replace the real anisotropic structure with 
isotropic equivalent, described by the equivalent dielectric 
constant £q, Or to use the “anisotropic material” option of the 
modern 3D simulators, introducing the pair of values &,,, and 
Eperp. The concept and design accuracy applying the first 
option was already described above. There exist only a few 
papers, where the second option has been directly applied, 
but with not so expressed benefits (see references in [A47]). 


Let’s compare both options. Several 3D models have been 
constructed for popular planar transmission lines on an 
anisotropic substrate, following the accepted in the RF 
design standard principles for ports, radiations boxes, etc. 
(see the example for MSL in Fig. 3.12a). We selected one 
typical Substrate 3.38 with a middle anisotropy, measured by 
the two-resonator method in the range 5-40 GHz with 
parameters &ar = 3.606; & erp = 3.345. The corresponding 
equivalent dielectric constant is &4 = 3.48 (design value), 
measured by the DPL method (catalogue value is 3.38). Fig. 
3.12b presents the simulated frequency £p -dependencies of 
a 50-Ohm MSL in two options of the used 3D simulator (in 
this case ANSYS HFSS): 1) isotropic equivalent substrate 


3.8 - - 
Substrate 3.38 











0 5 lO 15 20 25 30 35 40 
a f, GHz ©) 
Fig. 3.12 (a) 3D model of MSL on an anisotropic substrate with selected 
two rectangular volumes S (substrate) and A (air above) and two xOz planes 
Ps and Pa at the input port; (b) Comparison between simulated frequency 
dependencies of the effective &p and the corresponding equivalent &, 
dielectric constants of MSL on an anisotropic substrate using simulator 
options “isotropic or anisotropic material” 
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Fig. 3.13 (a) Measured and calculated £, values by simulated phase delay 
in several planar transmission lines and simulated normalized parallel and 
perpendicular E fields in these lines; (b) Distribution of normalized parallel 
and perpendicular electric fields in the united box volumes S+A (substrate 
and air box above) of 50-Ohms MSL (dashed curves — for united Ps + Pa 
planes on the radiation box wall — Fig. 3.12a) 










with & = S) = & = &, and 2) anisotropic substrate with uni- 
axial anisotropy, namely & = &) = Gari Ez = Eperp. The Eef 
values calculated by the second (anisotropic) option are 
visibly smaller, than the values for the isotropic case (the re- 
calculated equivalent dielectric constant is now 3.38 instead 
3.48 for MSL, 1.e. 2.8-% decrease). Fig. 3.13a presents these 
differences for other investigated planar lines: from 0.5 up to 
3 % decrease, which is not acceptable for a reliable RF 
design. In fact, the biggest differences appear for MSL and 
MSL-like lines, e.g. for coupled MSL (up to 4 % decrease, 
not shown). In the same time the course of the curves, 
representing the distribution of the normalized parallel and 
perpendicular electric fields, Epa Eror and Eperpl Eron in the 
whole structures fully correspond to the position of the 
measured &, values between &ar and &erp values of each 
given planar line — from Fig. 3.9c. Therefore, we can try to 
bind the ¢,, values directly to the field distribution. In [A46] 
we proposed a more accurate, field model for calculation of 
the €, values. For this purpose, we include in the MSL 
model (Fig. 3.12a) several new artificial volume and plane 
objects, where we can easily calculate the integrated 
complex E fields by the incorporated “field calculator” in the 
simulators. The aim is to determine both field ratios Epa Erot 
and Eperpl Erot 1n these volumes or planes and to calculate 
Eeq_MmsL DY the expression 
E 


par 


E x E par 
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The obtained results are quite informative. The calculated 
ratio for substrate only is Ena/Epery = 0.19/0.81 in the volume 
S. If we determine &  usz by (11) applying these ratios, the 
obtained dependencies are very close to the &g sim phase 
dependencies, obtained by the simulated phase delay (1.e. by 
the MSL S parameters) — see Fig. 3.12b. But if we now apply 
the E-field ratios in the whole structure, 1. e. Ena! Eperp = 
0.39/0.61 in the united volumes S+A, the calculated &, 
values are very close to the measured values & 9 meas (Or tO 
the values, obtained by the simulator for option “equivalent 
isotropic substrate”). This fact is very optimistic. It points the 
direct connection between the calculated actual ratio 
EparlEperp for each planar structure (MSL in the considered 
case) and the measured equivalent dielectric constant. 


The obtained positive results for the possibility to predict 
the actual values of the equivalent dielectric constant only on 
the base of the calculated near-field distribution of the 
parallel and perpendicular electric fields in MSL allow us to 
propose a reliable method for a direct determination of this 
important parameter without to perform any specialized 
time-consuming measurements. But to implement this 
method, we have to investigate the influence of the MSL 
parameters over the field distribution. Fig. 3.13b presents the 
most important dependences Epg/Eio: and Eperp/Eror on the 
substrate thickness, which can be used directly in expression 
(11). According to the expectations, the parallel component 
of the electric field plays a bigger role for thinner substrates 
and v.v. Therefore, we can average the dependencies and 
write an approximate expression 
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(11.1) 
which ensures the best accuracy for substrates with thickness 
20-30 mils. The other MSL and substrate parameters (e.g. 
metallization thickness, absolute dielectric constant, losses) 
don't play so important role. Now, if we calculate values 
Eeq_MSL_cae by (11.1), they practically coincide with the 
measured ones &g msz_iso IN the presented case in Fig. 3.12b. 


Caa (0.38 +0.02)x £ 
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IV. DIELECTRIC PROPERTIES AND ANISOTROPY OF 
REINFORCED, IMPREGNATED, WOVEN, KNITTED AND OTHER 
FIBRE MATERIALS 


The dielectric mixtures based on two or more dielectrics, 
one part of them performed as reinforcing filaments (fibres, 
threads, yarns) and another part — as appropriate fillers 
(ceramics, PTFE, epoxy resins, air, etc.), are very popular in 
the modern electronics — see Fig. 4.1. This is a wide group of 
useful artificial materials, which can be used as substrates for 
printed circuit boards (PCB), reinforcing skins and foam 
cores for multilayer antenna radomes, textile fabrics, used for 
wearable antennas, self-reinforcing plastics, etc. The most 
spread of them is the group of reinforced PCB substrates as 
commercial products (considered in IVA). They consist of 
woven or nonwoven reinforcing fibre-glass fabrics and 
different fillers and the resultant homogenized structures can 
be used in different parts of the microwave range. Nowadays 
some traditional textile fabrics also can act as substrates for 
the popular wearable antennas in the new communication 
body-area networks ([VB). Another variant of the reinforced 
epoxy-glass laminates is applicable for antenna radomes and 
other radio-transparent and mechanically stable constructions 
(IVC). Multilayer and gradient absorbers (IVD) are artificial 
materials, which could be added to this group. 


Of course, the question that must be answered when 
designing such composites is how to predict the effective/ 
equivalent dielectric properties of the whole system? (Note 
that similar approaches should be applied also to the 
mechanical and thermal properties.) We already considered 
in IIB several closed-form expressions for analytical 
determination of the equivalent complex dielectric constant 
of different mixtures (early concepts [9] and some new 
approaches for 2D [40, 41, 42] or 3D [43] woven fabrics, 
including the numerical techniques [44]). To obtain the 
required information for the resultant dielectric mixture, it is 
necessary to estimate the permittivity of the separate 
constituents (phases) and the used mixture topology. Finally, 
appropriate measurements of the resultant complex dielectric 
constants are also necessary to be performed by destructive 
or nondestructive methods (some of them considered in III). 





— 


B (a) courses —> (b) 





Fig. 4.1 Examples for some filament mixtures: (a) 2-D woven fabrics; (b) 
Jersey knitted fabrics; (c) non-woven fabrics; (d) 3D woven composites 
(legend: 1 — warp yarn; 2 — weft yearn; 3 — Z-yearn) 


The topology of the considered reinforced dielectric 
mixtures (shown in Fig. 4.1) suggests the existence of 
artificial anisotropy, caused by the spatial inhomogeneity in 
the filling. We already discussed the anisotropy origin (IC) 
and developed authorship methods for numerical (IID) and 
experimental (IIIB) determination of this parameter for such 
materials. We use as a measure for the degree of this 
anisotropy (separately for dielectric constant and dielectric 
loss tangent) the considered above parameters AA, AAtanés 
defined by expressions (9.1, 9.2), where &,, Sy and &, are 
the real dielectric constants, while tand,o,, tand,oy and 
tano,9, are the dielectric loss tangents (dissipation factors) 
along the axes Ox, Oy and Oz (as shown in Fig. 4.1d). The 
considered case refers to the biaxial anisotropy (for a pair of 
three different components along each axis). When the 
parallel components according to substrate surface coincide 
(i.e. &,=Hy=Epar ANd H-= Eerp; tandgo,= tanOzoy= tan Ôzpar and 
tanO,0,= taNOzye,p), the uni-axial anisotropy takes place. In 
the next four sections, we will consider the specific features 
and peculiarities of the considered types of dielectric 
mixtures and their anisotropy. 


A. Reinforced Commercial Substrates 


The reinforced substrates are key materials in modern 
electronics — the basis for the realization of new devices by 
microwave integrated circuits (MICs) and especially for the 
last SG communication standard. These materials consist of 
reinforcing fibre-glass fabrics and a variety of fillers — see 
Fig. 4.2. This mixture is a classic example for uni-axial 
anisotropic materials: usually, the glass fibres are high- 
permittivity supports, while the fillers — low-permittivity 
homogenizing media (with some exceptions, when the fillers 
are high-permittivity ceramic powders). Due to this 
combination of constituents the parallel dielectric constant 
(along the glass fibres' axes) is usually bigger than the 
perpendicular one, 1.€. Ear > Eperp (the behaviour of the textile 
fabrics with air filling is very similar; see IVB) 


The main role of the substrate is to support the metallic 
layout and surface-mounted components of the PCB schemes 
at low frequencies; however, in the microwave range, the 
substrate dielectric properties influence the electrodynamic 
behaviour of the whole printed structure. Thus, a reliable 
determination of the dielectric parameters of each 
commercial substrate by the manufacturers is very important 
for the consumers. The producers apply the reference IPC 
TM-650 2.5.5.5 test method [17] (based on clamped stripline 
resonator), which gives dielectric constant and dissipation 
factor along the Oz axis (Fig. 4.2a), denoted with the popular 
in the datasheets symbols Dk and Df and this pair of two 
values is usually fully enough for many applications at lower 


glass fibres 






0z _ filling 
microstrip 
Ox Lr conductor 
RAN pO Sol 
ARRRRRRRRRRRAAAO AAS 
SAANS, 
YY RNY = y- 
V Ren a : 
) SISOS = 
-o a (c) 


Pe ee a a a a a a 
x SS se ee ee ee ee ee 


ground conductor 


(a) (b) 


Fig. 4.2. Schematic view of the reinforced substrates: (a) substrate with 
uni-axial anisotropy [A21]; (b) top view (glass fibre fabrics and fillers); (c) 
side view of a microstrip line on a reinforced substrate 
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frequencies (for example, in planar schemes on the popular 
epoxy-glass laminate FR-4 [45], manufactured by many 
companies). However, at higher frequencies the anisotropy 
becomes important, especially for mmMICs, working in the 
mm-wavelength range, because the effective wavelength 
becomes compatible with the spatial irregularities in the 
reinforced substrate (well-illustrated in Fig. 4.2c). 

As we mentioned above (IIIB), we maybe were the first 
researchers (based on our experience with ferrites and 
ceramics), which detected the unwanted dielectric anisotropy 
of several popular commercial reinforced substrates in 2000- 
2002 by the developed authorship two-resonator method and 
published the first paper [A6]. For example, we measured for 
a commercial substrate Ro3203 (0.254-mm thick; 10 mils) 
considered as pure isotropic one (reference value 3.02 [46]): 
Epar ~3.238, Eperp ~3.036, equivalent dielectric constant &,~ 
3.111 in the X band and dielectric anisotropy AA,~6.4% 
(indeed the term “equivalent dielectric constant” was firstly 
introduced exactly in this paper). We presented also results 
for other popular substrates (e.g. for the practical isotropy of 
Ro3003, RT Duroid 6002; useful data for Ro4003, FR-4, 
etc.), but due to the used perturbation approximation for TM- 
mode resonator, some of the obtained results were not so 
accurate, especially for thick substrates. Then we developed 
more accurate full-wave analytical models and later on effec- 
tive numerical techniques of the measurement resonators and 
extraction procedures for determination of the parameters 
Dk/Df of uni-axial anisotropic single and multilayer samples 
by the two-resonator method ([A17, A21, A26, A31]). In a 
relatively short period (2004-2008) we succeeded to collect 
useful data for the actual anisotropy on many commercial 
substrates from different manufacturers and had the unique 
possibilities to compare and to publish some of them. In our 
practice, we met commercial substrates with different degree 
of anisotropy: 1) near-to-isotropic (AA, <2-3%; non-woven 
substrates); 2) substrates with middle anisotropy (AA, ~5- 
11%; the most typical case) and 3) relatively big anisotropy 
(AA, > 15-20%) [A32]. Today, we can distribute in each 
group many well-known commercial products on the base of 
our measurements, but this information is very sensitive for 
the most manufacturers and we try to use in our recent 
publications virtual denotations instead the actual substrate 
names when the products are still on market. The large PCB 
companies continue to support an opinion that the anisotropy 
of part of their products is “a discrediting property” for the 
consumers (although it is a completely natural property for 
such mixtures) and they usually don’t comment the measured 
anisotropy effects in details, nevertheless that the degree of 
this anisotropy directly depends on the applied technology. 
Until 2006 we didn’t manage to find any information about 
the anisotropy of commercial substrates published in the 
official specifications (application notes). However, in this 
period, 2002-2005, we found out several articles (not cited 
here), presented on the largest microwave conferences (e.g. 
IMS), where the authors disagree with part of the used 
reference Dk data. For example, some of them have assumed 
a value of ~3.52 for Dk of the very popular at that time 
substrate Ro4003 instead the catalogue value 3.38 (4.1-% 
increase!). They explained this decision with “better results 


obtained during the simulations of MSL filters on this 
substrate”. It was a strange explanation for us; RF designers, 
who typically draw the geometrical structures of the 
designed planar devices in the simulators with extremely 
large details, frivolously start to adjust the substrate Dk value 
in their projects until the simulation results coincide with the 
measurement ones (!). For information, in this period we 
measured the following Dk values for Ro4003 (0.508-mm 
thick; 20 mils): &ar ~ 3.66, & erp ~3.37 (close to catalogue 
value 3.38), equivalent dielectric constant exactly &,~3.52 
for 50-Ohms microstrip lines (MSL) in the Ku band; 1. e. a 
moderate dielectric anisotropy AA, ~ 8.3%. Other research 
groups also started to show their results for anisotropy of 
different substrates, mentioned in HC. 


Step by step, the large PCB manufacturers accepted the 
problems with substrate anisotropy and offered solutions to 
overcome them, but without to discuss in details the source 
of problems. The first company, according to us, which 
officially began to share information for measured dielectric 
anisotropy of its substrate, was Taconic Headquarters Ltd., 
Advanced Dielectric Division. In paper [47], the authors 
started to discuss the dielectric anisotropy of their material 
TLY-5A (measured by Bereskin’s method [48]). As early as 
2012, the next company that began to share results for the 
measured anisotropy of some of its products was Rogers 
Corp. Ltd. In the paper [49], the company researchers 
presented relatively detailed information for the anisotropy 
of some their substrates; even adduced a discussion for the 
obtained results, published earlier in our papers. They 
confirmed all our results, excepting data for one substrate 
with high-Dk. In 2013, Isola Group Ltd. also started to 
discuss and to solve the problems with substrate anisotropy, 
especially at higher frequencies (for substrates with 5G 
applications) and shared information for the equivalent 
dielectric constant [A40]. The concept for the equivalent 
dielectric constant, when the anisotropic substrate can be 
replaced with isotropic one, was more or less implemented in 
the datasheets. For example, Rogers Corp. (and then 
practically all other large PCB manufacturers: Isola Group, 
Taconic, ITEQ Corp., Shengyi Technology Co, etc.) started 
to give so-called design Dk value of each RF substrate 
(together with the reference technological Dk value obtained 
by IPC-TM-650 2.5.5.5 method) to ensure optimum simula- 
tion results. The large software developers also added the 
corresponding design Dk values in their material tutorials. 
Strictly speaking, the parameter design Dk was not well 
validated in the beginning, but the recent publications from 
the company researchers (e.g. [50]) discuss the problems of 
the substrate anisotropy and give unpublished data for many 
products. Nowadays, most of the companies present three Dk 
values (but not for all products): technological Dk (close to 
Eperp); design Dk (close to €, for MSL) and a value, obtained 
by SPDR method (close to &,,), which completely corres- 
pond to our concept for the reinforced substrates, developed 
since 2000 (see IIC). This information helps to solve the 
problems with substrate anisotropy; nowadays, the producers 
try to obtain “a full 3D dielectric portrait” of their products 
to have better control on technology, while RF engineers 
would like to perform a better 3D design of their devices. 
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Fig. 4.3. (a) Comparison between the frequency dependencies of 
equivalent dielectric constant &, for uncovered and covered MSL and CPW, 
compared with &, , and &¢, values of Substrate 3.38; (b) covered MSL line 


resonator (side view) for direct measurement of &,; (c) screened MSL line 
resonator (quasi-IPC method) for direct measurement of & erp [A44] 





Dielectric overlay 





Therefore, we can conclude that the development of 
measurement methods for determination of the important 
dielectric parameters of widespread commercial substrates, 
especially in the mm-wavelength range, 1s substantial and we 
continue to develop new methods (III). For example, in our 
paper [A44] we have presented new methods for direct 
measurements of the two important substrate parameters by 
direct determination of effective &: design Dk (by covering 
of MSL with enough thick overlay from the same material) 
and technological Dk (by an appropriate screening of MSL) 
— see Fig. 4.3. The measurements of close to parallel Dk by 
covering of CPW were also successful [A45, A58-A62]. The 
last new results are related to the temperature behaviour of 
substrate anisotropy, which are very interesting [A66]. 


Of course, our main advantage in the process for 
characterization of all accessible on the market commercial 
substrates is the possibility really to obtain a full 3D 
“portrait” of the dielectric properties of these important 
materials, measuring the anisotropy by our two-resonator 
method [A21, A26, A31] and equivalent parameters, when 
these materials have been considered as isotropic ones [A31, 
A42]. This information can help many users to perform 
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Fig. 4.4. Pairs of optimized measurement resonators with different 
diameters: (a) TE-mode R1 (30, 18.1, 15, 10 mm) and (b) TE-mode R2 (30, 
18.1, 10, 8 mm), which cover frequency range 5-39 GHz 
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Table 1. Measured dielectric parameters and anisotropy of several 
commercial substrates with practically equal or similar catalogue parameters 
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ee IPC TM 

Substrate h,mm_ parallel perpendicular equivalent ee , 6502.5.5.5 @ 

. EparltanSzpar — Eperp/tandeperp eq tan dz eq Mande, % 10 GHz 

Rogers Ro4003 0.510 3.67/0.0037 3.38/0.0028 3.53/0.0031 8.2/27.7 3.38/0.0027 
Arlon 25N 0.520 3.57/0.0041 3.37/0.0033 3.37/0.0033 5.8/21.6 3.38/0.0025 
Isola 680 0.525 3.71/0.0049 3.32/0.0042 3.32/0.0042 11.1/15.4 3.38/0.0030 
Neltec NH9338 0.520 4.02/0.0051 3.14/0.0025 3.51/0.0032 24.6/68.4 3.38/0.0025 
Rogers Ro3003 0.27 3.00/0.0012 2.97/0.0013 2.99/0.0013 1.0/-8.0 3.00/0.0013 
Rogers Ro3203 0.26 3.18/0.0027 2.96/0.0021 3.08/0.0025 7.2/25.0 3.02/0.0016 
Neltec NH9300 0.27 3.42/0.0038 2.82/0.0023 3.02/0.0023 19.2/49.2 3.00/0.0023 
Arlon DiClad880 0.254  2.32/0.0016 2.15/0.00093 2.24/0.0011 7.6/53.0 2.17/0.0009 





better and reliable design of different planar structures on 
these substrates, especially in the mm-wavelength range and 
the manufacturers to have better technology control of their 
products. Fig. 4.4 presents the photography of a set of eight 
TE- and TM-mode resonators, which “covers” the frequency 
interval 5-39 GHz by applying the fundamental modes in 
each resonator (and even up to 80 GHz by applying selected 
appropriate high-order modes). Using this measurement tool 
and its modifications, we managed to characterize more than 
70 different substrates, accessible on the world market; a 
small part of the obtained results (<10 %) have been 
published with the actual substrate names (e.g. as in [A31]) — 
as in Fig. 4.5 and Table 1. The other unpublished part 
composes our rich database for many substrates, including 
the last available products. We can characterize the separate 
substrate materials in several ways [A36, A40]. The first two 
ways are based on obtaining information for the anisotropy 
of commercial substrates versus the frequency (as in Fig. 
4.6a,b) and thickness (as in Fig. 4.6c,d or Fig. 4.55). The set 
of averaged values of parameters £ar, Eperp) Eeg aNd tan Og par, 
tan Oz perp, LANOzeg Zives valued information for each substrate 
and allows comparison between substrates from a specific 
group (like examples in Table 1), which is important for the 
producers and useful for the users. The last two ways for 
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Fig. 4.7. Statistical data for substrate NH 9338 (a) and its parameters’ 
inhomogeneity (b) [A13] 


Table 2. Measured standard deviations for parameters of large-size 
substrate sheets [A13] 


Substrate 


Ro4003 
NH9338 


SD Epar SD Eperp 


+0.5 
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32 


90 
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Fig. 4.8. Examples for measured parallel, perpendicular (&a;, &erp) and 
equivalent dielectric constants (€,) of different commercial substrates: (a) 
for different products of one company; (b) Comparison: case 1) low-Dk 
(~3.0) substrate with minimized anisotropy for small thicknesses; case 2) 
two low-Dk substrates with different anisotropy; cases 3, 4) family of 
substrates with Dk~3.38 from two different companies covering a wide 
thickness range 3-60 mil [A46]. Legend: ¢, measured by the DPL method; 
EeqCovMsL aNd Eyerpscrust. Measured by the direct method presented in [A44, 
A45] (see Fig. 4.3); &¢_caic calculated by the expression (11.1) [A46] 
obtaining of substrate parameters are difficult — statistical 
data, extracted from a large number (50-100) of samples (as 
in Fig. 4.7a) and information for dielectric parameters 
inhomogeneity in large-size substrate sheets (Fig. 4.7b). 
Such data, compared for a selected extract of several 
concurrent substrates, are decisive for the right choice of the 
larger users of these materials because they show the 
stability of the dielectric parameters and parameters of the 
constructed on their base devices — e.g. impedance Ze, 
propagation constant J and attenuation a@ (see Table 2). Fig. 
4.8 illustrates well how we can use this database for substrate 
parameters. A set of parameters and anisotropy of different 
30-mils substrates of one producer is presented in Fig. 4.8a 
with typical dielectric constants 3-3.8. Fig.4.8b (cases 1, 2) 
illustrates the attempts of another producer to minimize the 
anisotropy of low-Dk substrates with applicability in the 
mm-wavelength range; while Fig. 4.8b (cases 3, 4) gives a 
useful comparison between a family of products of two 
companies with similar Dk ~ 3.38 and different thicknesses. 


B. Textile Fabrics for Antenna Applications 


The textile fabrics show quite similar behaviour of the 
dielectric properties and anisotropy like the reinforced 
substrates due to the structures’ resemblance [51]. In our first 
publications on this topic [A47,A49], we presented measured 
by the two-resonator method dielectric anisotropy of several 
most popular natural and synthetic textile fabrics (see also 
selected results in Table 3). It turned out that the actual 
anisotropy AA, of textiles is small (1.5-5%) or moderate (5- 
9%) as the commercial microwave reinforced substrates 
(compare with data in Table 1). Only for more complex and 
multilayer artificial fabrics, the anisotropy AA, exceeds 10- 
12%. Then, in paper [A49] we developed 3D models for the 
numerical determination of dielectric anisotropy of woven 
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and knitted textile fabrics (Fig. 4.9) and presented useful 
results (part of them shown in Table 4). The information for 
the anisotropy of these artificial fabrics confirms the rule that 
the dielectric constant of these air-filled dielectric mixtures 
increases in the direction parallel to the axes of the fibres and 
decreases in the perpendicular direction (that why, different 
sample anisotropy AA, appears from +1-5% up to 12-18% — 
last column in Table 4) (for the promising flexible PDMS 


Table 3. Measured dielectric parameters and anisotropy of some classical 
textile fabrics (averaged values for frequency interval 0-36 GHz) [A62]. 
Insets at right: different disk textile samples 





























: 
> aiii 6 Epar Eperp | Anisotropy 
Textile fabric diii idr tand,pep | 44e/Adtianea% 
Epoxy-based =| 9 35 | 1.97/0.010 | 1.83/0.007 7.4/30 
waterproof fabric 
Waterproof fabric 
with breathability | 0.20 | 1.53/0.006 | 1.38/0.004 10.3/28 
GORE-TEX® [16] 
Weaved silk 0.19 | 1.60/0.028 | 1.54/0.016 3.8/57 
Weaved linen 0.65 | 1.65/0.043 | 1.58/0.044 4.3/-2.3 
Weaved hemp fishnet| 0.81 | 1.63/0.072 | 1.43/0.034 13.1/72 
Natural leather 0.84 | 2.47/0.055 | 2.44/0.054 1.2/1.8 
Weaved wool 2.10 | 1.28/0.026 | 1.21/0.015 5.6/54 
Jersey knitted wool | 5.50 | 1.40/0.024 | 1.26/0.021 10.5/13 
Denim 0.93 | 1.69/0.027 | 1.61/0.030 4.8/-11 
Cotton satin 5 0.25 | 1.58/0.019 | 1.45/0.013 8.6/38 
Jersey knitted cotton | 0.40 [ 1.56/0.055 | 1.50/0.044 He 9/22.2 





Fig. 4.9. Top and side views of 3D models of four artificial materials, 
made by straight or meander yearns with diameters 0.5 mm and distance 
between their centres 1 mm (the unit cell of samples is 1 x 1 x 1.5 mm; the 
unit cell has been repeated only in the Oxy plane) 


Table 4. Simulated dielectric parameters and anisotropy of several 
artificial textile structures 








i inii Epar Eperp! Anisotropy 
No Sample description /tan 5c par ian MA, /AAvanse,% 
1 Straight cylinders along to Ox O ai 1a 17.2/120 
direction (fish net) 
2 Meander cylinders along to 0x and 0Y 673/9.9022 1.750/0.0025  -4.5/-12.4 
directions (woven symmetrical net) 
s Meander cylinders along to 0x direction naa ; 5a 
and straight cylinder along to Ov directions baronet enamine PESAR 
4 Jersey knitted fabric 1.726/0.0026 1.530/0.0015 12.0/53.7 
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Fig. 4.10. Measured normalized resonance frequencies of length-(a) or 
width-bent (b) rectangular resonator 36x27 mm on substrates S1, S2, S3 and 
S4 (c) with measured anisotropy AA, ~ 0.4; 5.7, 13.7 and 11.3 % 
correspondingly [A62]. 





(Polydimethylsiloxane) polymer substrate, AA, ~ 6.4%; more 
details are given in in [A70, A64, A69]). These results give 
ideas on how to minimize the undesired anisotropy of these 
materials. In [A49] we investigated the behaviour of 
microstrip patch antennas and arrays on anisotropic textile 
substrates (frequency band, beamwidth, directivity and 
efficiency). In the recent paper [A62], we continued this 
research with bent microstrip rectangular resonators on 
selected substrates with different anisotropy. We found out 
that the resonance frequencies of the dominant mode in these 
bent structures (Fig. 4.10) depend not only on the curvature 
diameter Dc (obtained for pure isotropic substrates like S1) 
but also on the concrete anisotropy AA, (for samples S2, S3, 
S4) (or combination between Dc and A,), which is a comple- 
tely new result (the research is still in progress [A71]). 


C. Multilayer Antenna Radomes 


Antenna radomes are very important structures in the 
antenna design; they ensure the mechanical and weather 
stability of different open-air antennas and, at the same time, 
they should be radio-transparent in the operating frequency 
ranges. Therefore, the proper radome design separately from 
(or together with) the antenna design is an important task, 
because this structure adds losses and phase delay and 
influences the antenna beam-forming and directivity/gain. 
There exist two main types of radomes: thin (with thickness 
tr smaller than the wavelength 4) and resonance radomes 
(when tr ~ 4/2). The modern radomes are usually multilayer 
(sandwich-type) structures [52]. They consist of relatively 
thick middle core and two thin skin layers, which ensure the 
mechanical hardness. The cores are foams, wadding with air 
pores, honeycombs or similar low-Dk layers. On the contra- 
ry, the skins are reinforced layers: glass fabrics with resin 
epoxy filling, both with higher Dk. To perform reliable rado- 
me design, the dielectric properties of all radome layers have 
to be well determined. However, the considered structures of 
cores and hard skins show that they are more or less aniso- 
tropic materials. Exactly considering these needs to determi- 
ne the own anisotropy of radome layers (e.g. due to the 
destruction of antenna circular polarization at big elevations 
[A65]), we developed the authorship two-resonator method 
for multilayer samples, as it has been described in [A21]. 


Fig. 4.11 presents three examples of antenna radomes. The 
single-layer reinforced thin radome (Fig. 4.1la) could be 
considered as a reinforced substrate (IVA). The popular 
three-layer honeycomb radome consists of rare honeycomb 
cells (e.g. Kevlar® paper) and two different (top and bottom) 
skin layers — thin glass-fibre fabrics with resin filling — Fig. 
4.11b. All these samples have expressed anisotropy, 
illustrated in Fig. 4.12b. Based on these data we can 
construct different radome models — with one, three or six 


= 


(a) | oe b) e -AAS 

Fig. 4.11. Samples from radomes: (a) reinforced TWINTEX® fabrics 
woven with commingled E-glass and polypropylene filaments; (b) 
honeycomb radome: painted skin layer, Nomex® Kevlar-paper core and E- 


glass skin with glue residues [A19]; (c) Parabeam® 3D glass fabrics: 
reinforced 3D woven filaments with resin filling [A25] 
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Fig. 4.12. (a) Illustrative pictures of the used radome models with 1, 3 
and 6 layers; (b) dielectric anisotropy of separate layers and whole sample 
for each model; (c) IL-dependencies versus the incident angle @at 12.5 GHz 
for multi-layer honeycomb radome with total thickness 5.5 mm for normal 
(NP) and parallel (PP) polarization of the incident wave 
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Fig. 4.13. Dependencies of the insertion losses IL in Parabeam® 3D glass 
radome [A25]: (a) versus the incident angle at 12.5 GHz; (b) versus the 
frequency at several incident angles. The results are compared with Ansoft® 
HFSS simulations for 1-layer and 3-layer anisotropic radome models. The 
measured data are obtained by measurements of the received signal in DBS- 
TV steerable antenna array with and without covering radome in Ku-band. 
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Fig. 4.14. (at left) Half part of a large simulated structure: radome sheet 
between two open-end waveguides; (at right) Incident (a) and total E field in 
PML box with radome, illuminated by plane wave — normal NP (b, c, d) and 
parallel PP (e) polarization 
layers (Fig. 4.12a) (the last model includes two additional 
thin layers formed by the glue residues between the core and 
skins and a coating layer). We developed in our papers [A14, 
A19] an analytical model for determination of radome 
insertion (IL) and return (RL) losses of multilayer samples 
with anisotropy on the base of Paris’ model [53] for isotropic 


multilayer radomes. The comparison (Fig. 12c) between the 
numerical and measured IL fully confirms the validity of our 
theoretical radome model. 


The third example is very interesting from a measurement 
point of view — Fig. 4.11c. We proposed for the first time in 
[A21, A25] the commercial ParaBeam® 3D glass fabrics to 
be used as an antenna radome. The specificity here is the fact 
that the radome body is single structure: 3D glass fabric 
impregnated with epoxy resin and the skin layers and core 
are formed by a single process. Therefore, we cannot directly 
measure the core parameters. For this purpose we presented 
in [A21] a de-embedding procedure for the extraction of 
dielectric parameters of the middle layer: 1) characterization 
of parameters of the whole sample; 2) characterization of 
both skin layers and 3) extraction of the parameters of the 
inaccessible for direct measurements 3D glass core. The 
measured layer anisotropy is high: AA, ~27% for the top 
smooth skin layer; ~45.7 % for the rough bottom skin layer 
and —5.6 % for 3D glass core. With the help of the 
constructed radome model, we determined the IL in radome 
for an incident wave with normal polarization (Fig. 4.13) 
versus incident angle and frequency. Simulated and 
measured losses practically coincide. In [A25], for the first 
time we presented also well-developed 3D models for 
effective simulations for multi-layer anisotropic radomes — 
see illustrations in Fig. 4.14. The free-space model between 
two rectangular open-end waveguides gives applicable 
results but by very time-consuming simulations. On the 
contrary, stylized models with plane-wave source allow 
faster simulations at different incident angles with 
satisfactory accuracy; however, an extraction procedure has 
to be applied for the determination of the insertion losses and 
phases due to the multilayer radome [A25]. 


D. Multilayer and Gradient Absorbers 


Other important materials widely used in the microwave 
range are the microwave absorbers. Their main role in the 
structures, where they have been incorporated, is to suppress 
the parasitic interferences, cross-talk interactions and 
reflections. Recently, the microwave engineers started to 
design the whole microwave devices together with 
incorporated absorbers; therefore they need to know the 
actual EM parameters of these materials. Usually, the 
microwave absorbers have been offered as foams, rubber 
sheets, brushed coatings or thin films with absorbing 
inclusions (recently absorbing metamaterials). They consist 
of carbon, carbonyl iron, ferrites and other high-loss 
inclusions; therefore, they may have as pure dielectric 
properties (like the carbon-containing absorbers), as well as 
magnetic properties [54]. In [A10, A15, A16] we presented 
our attempts to determine the dielectric and magnetic 
parameters of thin microwave absorbers. Unfortunately in 
this period, the commercial absorbers’ producers presented 
in the corresponding catalogues only data for the absorbing 
abilities (through the sample in dB/mm), but not the actual 
dielectric and magnetic constants and loss tangents. In our 
review paper [A27] we summarized and compared the 
possible methods for characterization of microwave 
absorbers. There are no problems with pure dielectric absor- 
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Fig. 4.15. (a) Measured dielectric and magnetic constants (6, 4) of two 
commercial absorbers by perturbation method; (b) measured product & x 4, 
by three different methods 


Table 5. Extracted values of the dielectric and magnetic parameters of 30- 
um thick nano absorbers (AUT) mixed with protective lacquer L (nC — 
consists of nano-carbon particles of diameter ~3-4 nm; nF — nano-ferrite 


Fe30; particles of diameter ~20-30 nm; nFnC — a mixture nF: nC=1: 1) 
AUT +L Gpar/tandy par (7, GHZ)  Eperp Itan ó; perp (f, GHZ) Lar (f, GHz) 
7.65/1.6 (12.8); 6.82/5.4 4.20/0.22 (7.6) 4.04/0.15 

(22.1); 3.49/5.5 (32) (17.5); 3.56/0.22 (27.4) 
3.75/0.58 (12.8);3.73/4.1 2.05/1.93 (7.6); 2.54/0.55 


Lherp (f, GHz ) 
nC 


0.93 (15.7); 0.51 (24.4): 


nF 1.06 (20.8) 
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Fig. 4.16. (a) Carbon-gradient sample ECCOSORB® HR A15; (b) 
dependencies of the dielectric constant and dielectric loss tangent of four 
disk slices from 15-mm thick absorber 
bers; it is difficult to separate the magnetic and dielectric 
properties of magneto-dielectric absorbers. The separation is 
relatively easy [55] by measurements of the resonance 
parameters of exited TEj9, modes in rectangular resonators 
with small prism samples; first modes with odd p-values (1, 
3, 5) will ensure extraction of dielectric parameters &,, tano,, 
while first modes with even p-values (2, 4, 6) — magnetic 
parameters yw, tano, Fig. 4.15a presents measured 
equivalent dielectric ¢ and magnetic u, constants by 
perturbation methods, while Fig. 4.155 compares the product 
E X L, obtained by different measurement methods. As an 
example, we give the parameters of commercial rubber 
absorber FM8 (Collector Magma, Slovenia) in Ku band: 6, 
~18.0; uw, ~1.5, tand, ~0.2; tano, ~1.0. Thick absorber 
samples cannot be well characterized by resonance methods 
due to the big losses. However, we proposed in [A61, A67] a 
procedure for separate extraction of dielectric and magnetic 
parameters for thin nanoabsorbers by the two-resonator 
method using TE and TM modes with maximums of electric 
or magnetic fields in the place of the sample. Table 5 
presents results for nano-carbon and nano-ferrite absorbers 
(Ear, taNOzyar have been extracted by the help of modes 
TEomn With m = 1, 2, 3 and n = 1, 3; erp, tands perp — Modes 
TMomno with m = 1, 2, 3; far — modes TMo,,,, with m = 1, 2 
and n = 2; perp — modes TEon, with m= 1, 2, 3 and n = 2, 4. 
Dielectric parameters of carbon-containing absorbers 


strongly depend on the frequency, which corresponds with 
the results from other methods. 


The last example is connected with the characterization of 
an innovative type of foamed absorbers with gradient carbon 
distribution — see illustrations in Fig. 4.16a. We performed 
relatively thin slices from the thick materials t ~ 14 mm and 
managed to confirm this distribution as for the dielectric 
constant s. (d) = 1 + (d/0)’, as well as for the dielectric loss 
tangent tan ó;, (d) = 0 + (d/ Y (dependencies in Fig. 4.16b). 
Such absorbers have considerable decreased reflectivity, we 
proposed in [A50] this material to be used for decreasing of 
the radar cross-section (RSC) of unmanned vehicles. 


V. DIELECTRIC PROPERTIES AND ANISOTROPY OF 
CRYSTALLINE MATERIALS 


As we mentioned in HC, the crystalline anisotropy of 
single- or poly-crystalline materials (glasses, ceramics, 
artificial soft and low-temperature co-fired ceramics LTCC, 
liquid crystals, ferrites, semiconductors, etc.) is one of the 
oldest known types of well-expressed material anisotropy. 
Typically, these materials are homogeneous, but the 
anisotropy appears due to the existence of different 
crystallographic axes in their lattices and the fact that 
charges oscillate by different manner along these directions. 
Our first investigations were connected exactly with this type 
of anisotropy of ceramic, ferrite and semiconductor disks, 
cylinders, rings and rods by different microwave methods 
[A1-A3, A5]. Then we continued with the characterization of 
crystalline and crystalline-like materials by the two-resonator 
method [A23, A34, A53, A56] and applying combinations 
with other broadband methods [A59, A60]. Here we will 
present short summarized information for these difficult for 
characterization materials. 


A. Crystals, Ceramics, Glasses, Semiconductors 


The determination of dielectric properties of single- or 
poly-crystalline materials is a classical problem in the micro- 
wave frequency range [56]. There exist different methods for 
characterization of the dielectric properties of crystals; the 
most accurate are the resonance ones [20]. However, an 
additional circumstance for these materials is the presence of 
relatively strong biaxial anisotropy — different dielectric 
parameters along the different axes, e.g. E&x # &) # &, and 
tandx, # tand,y # tand,,. In two our papers [A53, A56] we 
especially applied two-resonator method for characterization 
of high-Dk materials like ceramics, glasses, ferrites, semi- 
conductors, etc. These papers aimed to show the ability of 
this resonance method to determine the crystalline anisotropy 
with acceptable accuracy and how this information can be 
used for crystalline sample characterization. 


Table 6 (part 1) presents dielectric parameters and 
calculated anisotropy of different ceramics, performed as 
disk samples, which fit the measurement resonators’ 
diameters. They are low-loss and relatively high- materials, 
suitable for compact microwave integrated circuits MIC’s, 
but the relatively big difference between the absolute values 
of parallel and perpendicular dielectric constants is a serious 
design problem. Usually &par > erp; AA, ~ 9-15% (only for 


19 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


the popular polycrystalline Alumina substrate we measure 
Epar < Eperp; AA; ~ —7%). The attractive artificial soft ceramics 
(3M@®Epsilam10, Rogers® Ro3010, TMM 101), plastic 
substrates with high-e filling, show even bigger anisotropy 
(7-23 %), excluding the relatively isotropic RT Duroid® 
6010, AA ~4.4%). The new class of low-temperature co-fired 
ceramics (LTCC), applicable in the monolithic MIC’s and 
promising for the mm-wavelength range, also have 
noticeable anisotropy (AA, ~ 13%). In general, our results for 
high-g ceramics are accurate, but difficult to be obtained — 
disk samples with fixed large diameters should be prepared, 
which is technologically inefficient process. 


That’s why, we proposed in [A53, A56] another variant of 
the two-resonator method for easier determination of sample 
anisotropy: to use separately either TE- or TM-mode resona- 
tors (R1 or R2; IIIB), but now for smaller disk or prismatic 
samples with different orientations parallel or perpendicular 
to the E fields — see the illustrations in Fig. 5.1a,b,c. 

Table 6. [A56] (Part 1) Dielectric parameters and uni-axial anisotropy of 
ceramic samples, measured by two-resonator method (R1 & R2); (Part 2) 


Dielectric parameters and anisotropy of three disk samples measured by TE- 
mode resonator R1 in position D1 and D3 (resonator axis lies along 0z) 















No Sample; reference values ké] Anisotropy 
Part | Epar (tande par Eperp! tande perp AAs AAtande,% 
l Alumina (Al-0:}) 9.8-10.7 9.65/0.0003 10.35/0.0004 -7.0/-29 
2 Polycore (Al,O; + 0.3% MgO) 10.044/0.0002 9.21/ 0.0003 8.7/-40 
3 Sitall (glass ceramic) 8.19/0.0042 7.16/ 0.0038 13.4/10 
4 3M ®Epsilam 10 9.8/0.0020 11.64/0.0022 9.25/0.0045 22.9/-69 
5 Rogers ®TMM10i 9.9/0.0020 11.04/0.0019 10.35/0.0035 6.5/-59 
6 Rogers ®Ro3010 10.2/0.0035 11.74/0.0025 10.13/0.0038 14.7/-41 
7 RT Duroid® 6010 10.2/0.0023 10.71 10.252 4.4 
8 ACX® LTCC 7.5/0.003 7.60/0.007 6.68/0.0075 12.9/-6.9 
9 YIG Garnet 14.8; Ms 1.7 kA/cm 15.83/0.00020 13.11/0.00027 18.8/-30 
Part 2 Com 
la Alumina (A1,0;) 9.8-10.7 9.78 / 0.00080 10.20 / 0.00075 -4.2/6.5 
8a ACX® LTCC 7.5/0.003 7.66 / 0.0069 6.20 / 0.0048 21/36 
10 MgTiO; 14.61/0.00028 15.10 / 0.00018 -3.3/44 
Oz Oy 
L? f $ Pos AR 
fh it teh 1 s y 
Pos. P1 % R2 j: peii 
R i ee 
Pos. P2 Pos. P3 me 
(a) Oz 
Pos. D1 





Pos. D3 (c) 


Fig. 5.1. Possible sample positions in resonators for prisms (a) and disks 
(b); (c) photo of samples in TM-mode resonator R2 in positions D1 and D3; 
(d, e) illustration of the E-field vector distribution in the measurement 
resonators R2 (TMo;9 mode) and R1 (TEo;; mode) 


Table 7. Dielectric constants and anisotropy of several prismatic samples 
measured by TM-mode resonator R2 in positions Pl, P2 and P3 (resonator 
axis lies along with 0z) 


No Sample description on rI Cate] Aan% Ag% 
Exx Ey Ez: 

1 Quartz (single crystal, c-axis =0z) 3.74 3.84 4.27 -13.2 -10.6 
2 Quartz (fused) 4.32 4.31 4.325 -0.12 -0.35 
3 Silica glass 5.15 5.23 4.30 18.0 19.5 
4 Mica (multi-layer silicate) 4.87 4.70 3.00 47 44 

5 Optical glass LiNbO; 32.25 29.3 22.8 34 25 

6 Crystalline Si for wafer 10.75 10.15 8.125 28 22 

7 Semi-isolated GaAs wafer 11.30 11.10 7.90 36 34 





However, the analytical model of the measurement 
resonators, described in [A21, A26] cannot be used; suitable 
3D models for numerical simulations have to be applied now 
(examples have been given in Fig. 5.1d,e). Table 6 (part 2) 
presents again the dielectric parameters and anisotropy of 
some of the disk samples in part 1, but now measured by TE- 
mode resonator R1 in positions D1 and D3 — the coincidence 
is good. The variant with prismatic samples allows us to 
determine the actual bi-anisotropic parameters of different 
high- samples. Table 7 presents the measured values £x, Ey, 
and &, in positions P3, P2 and PI of several glasses and 
semiconductor wafers and the corresponding anisotropy 
AA xy. The new information from these data is the expected 
fact that the single crystals show bigger anisotropy than the 
anisotropy of the corresponding poly-crystalline sample. For 
example, fused and single-crystal Quartz samples have equal 
values for the perpendicular permittivity, £, ~ 4.27-4.32 (the 
c-axis coincides with the resonator axis 0z), while the values 
for the parallel permittivity &x, &y are different: ~3.74-3.84 
in the single crystal along the a-axes; ~4.31-4.32 in the fused 
Quartz. Relatively big anisotropy we observe for the crystal- 
line Si and semi-isolated GaAs (last two rows in Table 7). 


B. Microwave Ferrites 


The microwave ferrites, e.g. spinel, garnet, and hexaferrite 
systems as thin or thick films, powders, bulk samples, and 
nowadays different multiferroic materials, are key 
components in the systems that send, receive, and manipulate 
EM signals across a wide frequency range; from VHF up to 
mm-wave bands [57] and they can ensure non-reciprocal 
behaviour. Microwave ferrites are very similar to the 
microwave ceramics (high dielectric constant & and low 
dielectric losses), but with this difference that they 
additionally possess magnetic properties — small magnetic 
constant (yem <1) and bigger magnetic losses, when no 
external dc magnetic field Hp has been applied (we will not 
consider the tensor magnetic properties in non-zero field Ho). 
In our early papers [A2, A3, A5] we developed methods for 





* (O ee (d) 
Fig. 5.2. Measurement resonators and sample position: (a) Resonator R1 
with TEomn modes (m = 1, 2, 3, n = 1,2,3,4,5,6); (b) Resonator R2 with 
TMomo modes (m=1,2,3); (c, d) Photo of Rl and R2. Legend: Ih, lv — 
samples in horizontal and vertical position; 2 — resonator; 3 — support foam 





Fig. 5.3. E- and H-field distributions in resonators with samples in 
position lh and lv: (a) R1 with excited TEo;; mode; (b) R2 with excited 
TMo mode; (c) R1 with excited TEo12 mode 
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determination of the dielectric constant & and saturation 
magnetization M, of ferrite disks, prisms and cylinders, but 
calculated the demagnetized value Lem by the known 
expression Ladem =(1/3)[1+2(1-(m™M Jf"), where 7 = 3.5186 
MHz/(A/cm) is the gyro-magnetic ratio. Then in [A23], we 
began to measure the dielectric constant anisotropy of 
ferrites (Table 6, row 9). In our recent investigations of the 
Z-type hexaferrite Sr3Co2Fe.40,; we managed to measure for 
the first time the full set of dielectric and magnetic 
parameters and the corresponding dielectric and magnetic 
anisotropy (for Hop= 0) above 3 GHz [A60]. In this case, we 
apply the following measurement procedure. The disk ferrite 
samples have been measured by the pair of resonators that 
support either TEom, modes Gin R1) or TMomo modes (in R2) 
with samples placed in horizontal (1h) and vertical (lv) 
position in each resonator — see Fig. 5.2. The first several TE 
modes in R1 with m=1,2,3 and odd n=1,3,5 ensure extracting 
of parallel & ar, tanôzpar GN lh position) and perpendicular 
Eperp, tANOzyerp (Lv) dielectric parameters because the sample 
falls in the E-field maximums (Fig.4.19a), while the first TM 
modes with m =1, 2, 3 — extracting Of & erp, tan Ôzperp (1h) and 
Epars tANOs par (1v) (Fig. 5.3b). On the contrary, the first 
several TE modes with even n =2,4,6 in RI has minimums of 
the E fields and maximums of H fields in the place of sample 
and allow us to extract with enough accuracy the magnetic 
parameters: in parallel direction Jp, tanOypyar (IV) and in 
perpendicular direction per», tanOypery (lh) (Fig. 5.3c). In 
fact, only TE modes with even n index have pure H-field 
maximum and pure E-field minimums in the place of the 
sample (shown in Fig. 5.4; excluding very high samples, Fig. 
5.4b); the TE modes with odd n index and TM modes have 
mixed E and H fields in the place of the sample (illustrated in 
Fig. 5.3a,b). These fields exist in mutually perpendicular 
directions: modes, which allow 6&,,, extraction, have been 
influenced by 4perp and v.v. That’s why, the extraction 
procedure has been proposed as follows (the set of extracted 
values have been presented in Table 8). First of all, using TE 
modes with even n index and TM modes we determine 
intermediate values of the perpendicular and parallel dielec- 
tric parameters (Table 8: marked in grey) for w=1 and 


F-tield + Il-field 
N | 





Fig. 5.4. E- and H-field distributions: a) R1 with TE;;2 [A23]; b) R1 with 
TEoi2 mode; sample in position 1v: c) R1 with TEo2 mode; sample in 
position 1h 


Table 8. Extracted dielectric and magnetic parameters of hexaferrite 
Sr3Co2Fe2401; by the two-resonator method using different TE and TM 
modes for horizontally and vertically placed samples 


Dielectric parameters: 
e/tand, 
Par: 15.02/0.076 
Perp: 9.03/0.25 
Perp: 9.92/0.445 
Par: 14.76/0.026 


Magnetic parameters: 
ptand, 


Mode (resonance Sample 
frequency, Q factor) orientation 
TEo1; (12.4008/3 1.86) hor. 

(13.0374/85.08) vert 
TMo: (7.5620/31.86) hor. 
(6.9632/194.5) vert. 
TEo12(15.6616/550.6) hor. 
(15.6681/387.2) vert. 
TEo14(23.3829/963.1) hor. 
(23.3907/511.0) vert. 


Perp.: 0.92/0.70 

Par.: 0.70/0.32 

Perp.: 1.03/0.12 

| ) Par.: 0.73/0.40 
Par.: 14.7/0.08 | Perp.: 9.30/0.15 Par.: 0.75/0.20| Perp.: 0.98/0.10 
Dielectric anisotropy: 45/-61 % Magnetic anisotropy: -27/67 % 





Averaged values (5-20 GHz) 


tand, = 0. Using TE modes with odd n index we determine 
Llpar, taN Ôn par (1V) (for intermediate values & erp, tan Ôsperp) and 
Llperp) t€NOypyerp (1h) (for intermediate values &,,, tang par). 
Then, with the obtained intermediate magnetic parameters 
we extract the corresponding final dielectric parameters and 
finally, we apply the same procedure for the determination of 
corresponding final magnetic parameters (Table 8: marked in 
black). Of course, there exists some uncertainty for the final 
obtained € and zu parameters due to the differences between 
the resonance frequencies of the neighbour TE and TM 
modes. However, we apply additionally in [A60] a wideband 
method for final verifying of the frequency dependencies of 
ferrite material parameters (this procedure will be described 
later on). On the base of applied measurement methods we 
obtained the following averaged parameters for the 
considered Z-type hexaferrite Sr3Co.Fe24Ox4, in the frequency 
range 5-15 GHz: &4, ~ 14.7; tandgpar ~ 0.08; & erp ~ 9.3; 
tanó-perp ~ 0.15; par ~ 0.75; tanóppar ~ 0:20; Lperp ~ 0.98; 
tan Ô,perp ~ 0.10 (Table 8; last row). The measured anisotropy 
is relatively large, as for the dielectric parameters, as well as 
for the magnetic ones: AA, ~ 45 %; AA,, ~ -27 %. 


VI. PROPERTIES AND ANISOTROPY OF 3D-PRINTED 
DIELECTRICS 


The opportunity for 3D printing of different dielectric, 
metallic or mixed artificial materials provokes nowadays a 
lot of new applications of these materials in the microwave 
devices, including antennas. The 3D printing as a modern 
additive technology gives also the possibility to mix two or 
more dielectrics in one mixture with well-designed 
properties for microwave applications similar like the known 
reinforced substrates, textile fabrics, antenna radome 
composites, foams, absorbers, etc. Thus, the 3D printing is 
practically an ideal technology for producing of a variety of 
anisotropic metamaterials [14], but also isotropic or 
anisotropic filling materials with common purposes, e.g. 
foam-like and substrate integrated non-radiating dielectrics, 
and also metallized 3D printed devices. 


A. 3D Printed Foam-Like Isotropic Dielectrics 


We began to investigate in [A47] the 3D printed structures 
regarding possible artificial anisotropy of the resultant mate- 
rials, and established some preliminary conditions for obtain- 
ning artificial isotropy — see Fig. 6.1. The conclusion is that 
it is easier to produce anisotropic materials by 3D printing 
than artificially isotropic materials. For example, designers 
of Luneburg lens antennas by 3D printing [58] have to 
produce relatively large in size complex artificial dielectrics 
with a specific spatial distribution of the dielectric constant 
but have to ensure also a practical isotropy/ homogenization 
(Fig. 6.1, the unit cell 1). The main condition for obtaining of 
an isotropic distribution of the dielectric constant is the 
applying of fully symmetrical unit cells along all three 
directions, which usually is in a contradiction with reliable 
mechanical properties of the whole printed body. 


We discussed in our paper [A54] the possible 3D models, 
limitations and practical realization of artificial isotropic 
foam-like dielectrics by 3D printing, applying different 
suitable unit cells. The design strategy for the control of 
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Fig. 6.1. Calculated dielectric anisotropy (a) of several 3D printed 


dielectrics with different unit cells (b) versus the ratio between the dielectric 
constants of inclusions and hosting substrate 





dielectric constant anisotropy. % 


absolute values of the isotropic equivalent dielectric constant 
of these dielectrics has been demonstrated in Fig. 6.2 and 
Fig. 6.3. First of all, a suitable symmetrical unit cell can be 
selected with inclusions like cubes (Fig. 6.2a and [58]), 
spheres, prism, rhombs, cylinders or their combinations, 
supported by thin prismatic rods in all three directions. Then 
a bigger artificial substrate with appropriate dimensions has 
to be constructed by a repetition of the selected unit cell — 
Fig. 6.22b and Fig. 6.3. Finally, applying the procedure 
described in IID, we can calculate the equivalent (resultant) 
dielectric constant and loss tangent and compare with the 
measured values (Fig. 6.2c). The comparison for cube unit 
cell between our results and these in [58] 1s very good, at the 
comparison with measurement results for prism unit cell. 


This technology is very effective and accurate for 
construction of 3D printed homogenous dielectrics with the 
desired dielectric constant. It is ideal now for fast and accu- 
rate 3D printing on different lenses as antennas — [58, 59]. 
We developed 2009 in [A29] accurate step models for an 
approximation of the dielectric constant gradient distribution 
in the Luneburg lens antennas and introduced effective 3D 
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Fig. 6.2. Model of 3D printed sample: a) cubic unit cell 5x5 mm; b) 


whole substrate 15x20x5 mm; c) calculated and measured resultant 
dielectric constants by additive mixing technologies 
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Fig. 6.3. Artificial isotropic samples constructed by unit cells with: a) 
cube; b) sphere; c) prism; d) cylinder 








Ew 
Fig. 6.4. 3D printed samples by Formlabs2® stereolithography (SLA) 3D 


printer applying dental model photopolymer resin FLDMBEOl1 (see the 
measurement results in Table 9) 


Table 9. Measured dielectric parameters and anisotropy of 3D printed 
samples from Fig. 6.4 


Sample description thickness, mm Epar Eperp Anisotropy AA, 
1. Pure resin FLDMBE01 2.04 2.980 3.00 -0.67 
2. Cubes 2x2 mm 7.70 1.332 1.330 0.15 
3. Prisms 1x1 mm along to Oxy 9.30 1.497 1.450 3.19 
4. Prisms 1x1 mm along to 0z 9.33 1.320 1.445 -9.04 
5. Prisms 1x1 mm along to Oxyz 9.22 1.407 1.395 0.86 


models for simulations of several variants of these promising 
antennas. Now we realize that these early models are fully 
applicable in the 3D printing of lens antennas and we have 
renewed our research on this topic. Fig. 6.4 presents several 
new 3D printed samples with controlled dielectric isotropy or 
anisotropy. They have been produced by Formlabs2 stereo- 
lithography (SLA) 3D printer applying photopolymer resin 
FLDMBEO1, which ensures final size accuracy ~30 um. 


The measurements of the dielectric parameters of these 
samples by the two-resonator method allow us to prove the 
concept for the symmetrical unit cells and foam-like 
dielectric behaviour. Fig. 6.5 illustrates well this concept. We 
constructed basic 3D net with square prisms (side | mm) and 
added new prisms with different sides (2.5 and 4 mm) and 
tuned length, orientated in different directions. The result is 
an increase of the absolute anisotropy AA,, which depends on 
the degree of filling, nevertheless, if its sign is positive (for 
horizontal orientation) or negative (for perpendicular 
orientation of the prisms). Only for pure cubic unit cell, the 
anisotropy keeps its small absolute value (the presented 
research is still in progress). 


B. 3D Printed Dielectrics with Non-Symmetrical Unit Cells 


The 3D printers generally save time and material by 
making the interior of the printed object not completely solid 
(with low infill density). The hexagonal and rectilinear (grid) 
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Fig. 6.5. Dielectric constant anisotropy of 3D printed nets made by 
prisms with different sides and lengths (both normalized to the unit-cell 
length 5 mm) 
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infill 1s the most common because the samples are strength 
and light. However, due to the poor unit-cell symmetry, the 
resultant anisotropy of samples with hexagon infill is relati- 
vely large for electro-dynamic purposes. Fig. 6.6a presents 
the dependence of uni-axial anisotropy AA, on the outer 
hexagon side length a (normalized to the whole sample 
height A). The sign is negative for the chosen absolute height 
h= 2 mm; an absolute maximum appears, when the cell 
length is close to the height, i.e. a ~ h. In the case of mixed 
hexagonal-triangle filling, this maximum has been shifted for 
higher a/h. The anisotropy depends also on the width w of 
the hexagonal walls — see Fig. 6.6b; there also exists an 
absolute maximum for AA, (the sign is negative for the fixed 
height). However, the most interesting is the AA, dependence 
on the sample height h — see the corresponding curve in Fig. 
6.6b for the normalized height to the hexagon length a. We 
can see that the anisotropy has as positive, as well as a nega- 
tive sign. In our case, when the height h/ is about 20 % from 
the cell length a, the anisotropy is AA, ~ 0. This is an impor- 
tant fact — the equivalent dielectric constant of the considered 
3D printed dielectric depends on the sample height! It shows 
also how difficult 1s to achieve isotropy of such materials on 
the base of hexagonal (or other non-symmetrical) filling. 
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Fig. 6.6. Uniaxial anisotropy in 3D printed samples with hexagonal infill 
versus (a) normalized outer side length a (a/h); (b) normalized hexagon 
width w and height h (w/a; h/a) 





Fig. 6.7. 3D printed samples by fused deposition modelling with PLA 
filament applying low-cost 3D printer MakerGear® (a). Samples (b): pure 
PLA wire 1, samples 2-9 with infill 60-10 % and different height; dense 
(~100 %) 3D printed photopolymer 10 (see measurement results for samples 
1-10 in Table 10) 


Table 10. Measured dielectric parameters and anisotropy of 3D printed PLA 
samples from Fig. 6.7 


Samples h,mm Epar [tan dz par Eperp! tan de perp previa! 

1. PLA wire d=2 2.750/0.0080 2.750/0.0080 0/0 

2. infill 20 % 1.85 1.310/0.00260 1.273/0.00146 2.9/56 
3. 30% 1.69 1.368/0.00287 1.328/0.00161 3.0/57 
4.40% 1.88 1.555/0.00378 1.504/0.00229 3.3/49 
5.50% 1.42 1.690/0.00493 1.620/0.00281 4.2/55 
6. 60 % 1.68 1.810/0.00524 1.710/0.00310 5.7/178 
7.10% 6.70 1.177/0.00158 1.216/0.00177 -3.3/-11.3 
8.10% 1.30 1.111/0.00128 1.134/0.00093 -2.1/31 
9.10% 0.96 1.035/0.00095 1.071/0.00052 -3.4/58 
10. photopolymer 100 % 1.69 2.961/0.0049 2.960/0.0053 0.03/-7.8 





We were able to prove this effect experimentally. Table 10 
presents the measured dielectric parameters and anisotropy 
of several home-made samples from polylactic acid (PLA) 
polymer, applying a very low-cost 3D printer. 


The photos of some produced and measured samples have 
been given in Fig. 6.7. Nevertheless the non-perfectly 
produced samples, we can see the tendency: the anisotropy 
AA, increases with infill increase (3-6%; stronger for AAtané¢). 
A negative anisotropy has been observed for samples with 
predominant vertically deposited walls. A dense (~100-%) 
3D printed photopolymer sample shows. well-expressed 
homogenization, which is an expectable fact. All these 
investigations are still in progress. 


C. Anisotropy of Substrate Integrated Nonradiative 


Dielectrics (SINRD) 


A new class of substrate integrated nonradiative dielectrics 
(SINRD) has been proposed in [60]. They are ordinary 
substrates with a drilled pattern of air holes. The idea is to 
decrease the equivalent dielectric constant and dielectric loss 
tangent and to increase the applicability of these artificial 
materials in the mm-wavelength range. A similar approach 
can be used also for easy construction of Luneburg lens 
antennas, made by slices with incorporated air holes [61]. 
However, in these cases, we again detect the already pointed 
reasons for appearing of uni-axial anisotropy — directed 
inclusions (air cylinders in this case), placed perpendicularly 
to the isotropic substrate surface. We performed special 
research on the influence of the incorporated air holes of 
diameter d at distance L on the dielectric constant and 
anisotropy of isotropic polycarbonate (PC) substrate. As in 
the case of 3D printed materials with non-symmetrical (e.g. 
hexagonal) filling, we found out that the dielectric anisotropy 
of SINRD strongly depends on the air inclusions’ dimensions 
and substrate height h — see Fig. 6.8. The anisotropy AA, 
could be as with positive, as well as with negative sign for 
small L/d values. 
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Fig. 6.8. Dielectric constant anisotropy due to the air-filled cylinders in 


isotropic PC substrate (SINRD) 


Table 11. Measured isotropic PC and anisotropic Arlon 25N samples 
without/with air via holes 





Samples + ) h, mm Epar Itan; par  Eperp! tAn Ô; perp ag 
Polycarbonate (PC) mr = (0.50 2.765/0.0050  2.764/0.0049  0.04/2.02 
PC + air holes (d ~ 1.7) a 0.50  2.474/0.00428 2.421/0.00409 2.17/4.54 
Arlon 25N E 0.513  3.645/0.00401  3.309/0.00287  9.67/33.14 
Arlon 25N + air holes (d ~ 1.5) 725% 0.513  3.413/0.00365 3.113/0.00267 9.21/31.01 





We tested experimentally by the two-resonator method the 
influence of perpendicular air via holes in isotropic PC 
sample and anisotropic commercial substrate Arlon 25N 
(results in Table 11). The anisotropy AA, increases for 
isotropic PC with air holes to 2.17 %, while the anisotropy in 
the anisotropic substrate slightly decreases due to the 
incorporated isotropic (air) inclusions. 


D. Anisotropy of Fresh Plant Tissues with Cell Structure 


The dielectric anisotropy of the fresh plant tissues with 
cell structure can play an important role in the determination 
of the so-called fresh weight (FW). The nondestructive FW 
determination of young plant shoots, leaves, seeds, blossoms 
etc. is necessary for the diurnal regulation of plant growth 
and ensures valued information for the growing processes 
[62, 63]. In [A48] we managed to show that the dielectric 
anisotropy of these tissues is very large as in a fresh state, as 
well as in a dry state, probably due to the specific hexagonal 
cell structure (as in the case of 3D printed samples). Walnut 
leaves have been selected for the anisotropy test. We made 
disk samples from the leaves and put them in vertical or hori- 
zontal positions in resonators (Fig. 6.9a,b), which support 
simultaneously TE ,;; and TMo;o modes (Fig. 6.9c). Due to 
the different water content, the resonance shifts depends on 
the actual sample FW and can be used for building of a 
calibration curve “frequency shift Af, MHz — fresh weight 
FW, mg” by staking together a different number of equal 
disk samples, smoothly increasing FW. However, we found 
out that the calibration curves differ for samples orientated 
vertically and horizontally; dependencies are given in Fig. 
6.10. The reason definitely is the large leaves anisotropy. We 
determined the dielectric parameters of fresh and dry leaves 


— Table 12; both dielectric constants and anisotropy AA, 


decreases with the sample drying, remaining large — from 92 
up to 35%. Thus, if we apply one of the obtained calibrations 
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Fig. 6.9. (a) Walnut leaves and disk samples under test; (b) 3D models 
plant tissue of disk samples from leaves, placed horizontally and vertically 
on foam support in measuring cylinder resonator (c) with pair of TMoi9 and 
TEoıı modes (D = 150 mm, H = 172 mm) 


Table 12. Data for measured parallel and perpendicular dielectric 
parameters and anisotropy of samples from fresh and dry walnut leaves 











Sample description Epar [tAn z par Eperpl tan Sz perp ak 
Immediately after the pick-up of the leaves 24.8/0.315 9.15/0.235 92/29 
1.5 hours after the pick-up of the leaves 16.1/0.305 7.65/0.225 71/30 
Dry leaves 2.0/0.0525 1.4/0.0275 35/63 
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Fig. 6.10. Calibration curves for horizontal and vertical sample 
orientation: frequency shift v/s FW in mg for TMoj0-mode and positions on 
the plot for some unknown samples (single straight leaf, once or twice 
folded leaves, I -shaped or crashed walnut leaves) 


Straight leaf T-shaped leaf 


curves for FW determination of plant with different leaves 
orientations (as in Fig. 6.1le), we will get wrong results. To 
prove our concept, we constructed 3D models of a fictive 
plant organism with three leaves (with stem) in different 
orientations, which have noticeable dielectric anisotropy 
(parallel/normal dielectric constant 35/5) — Fig. 6.1 la-d. The 
obtained frequency shifts in MHz and recalculated relative 
FW are given in the insets. We can see that in both opposite 
cases (a, b) the shifts differ 19 times (!) for TMo 9. For leaves 
in an inclined position (c) or a mixed position (d) the shift is 
practically the mean value from both opposite shifts. For 
example, if we measure crashed or I’-shaped (half of the leaf 
is placed vertically and the other half — horizontally) walnut 
leaves, the point for FW fall exactly in the middle zone 
between the pair of reference curves in Fig. 6.10. This fact 
gives a possible solution for the right prediction of the FW of 
the whole plant organism by nondestructive resonance 
method — to introduce a correction coefficient K depending 
on the predominant leaves' orientation in the real plant. Let 
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Fig. 6.11.3D models of a fictive plant organism with 3 leaves placed: (a) 
vertically to the resonator axis; (b) horizontally; (c) with slope 45 deg; (d) 
mixed position (data in boxes: frequency shift and calculated FW relative to 
the case (b); (e) examples for plants with different leaves orientation 
the averaged relative portion of the horizontally placed 
leaves in the whole plant organism is x (determined by 
optical methods). If the frequency shifts for pure horizontal 
and pure vertical orientations of the plant leaves are Af, and 
Af,, the correction coefficient could be K = (1-x)Af, + xAf, 
and the actual fresh weight will be FW = K.FW,, where FW,, 
is the fresh weight, which corresponds to pure horizontal 
orientation. For example, a mean reference curve is given in 
Fig. 7 for x = 0.5 (50 %). Thus, for FW = 500 mg of the I- 
shaped sample (x = 0.5), the calculated coefficient is K ~4.61 
(the actual shift is 4.72, the measurement uncertainty is less 
than 2.4 %; instead ~600 % without taking the anisotropy 
into account). Therefore, the correction procedure allows 
enough accurate determination of FW taking into account of 
the actual plant tissue anisotropy. 


VII. DIELECTRIC PROPERTIES AND ANISOTROPY OF 
MICROWAVE METAMATERIALS 


The engineered microwave metamaterials and photonic 
band-gap crystals can be designed with controllable 
dielectric/magnetic constants along the different directions. 
In this case, the artificial anisotropy is a fully desired 
property, which ensures unusual characteristics of some 
anisotropic metamaterial devices [7, 15] like invisibility 
cloaks, EM concentrators, EM-wave converters, etc. We will 
concentrate our attention on specific issues connected with 
experimental dielectric properties’ characterization of these 
materials. As we already pointed in III, there exist a lot of 
measurement methods for material characterization, but only 
a few of them can determine with acceptable accuracy the 
artificial anisotropy. In this section, we will share our 
experience for anisotropy characterization of metamaterials. 


A. Metamaterials with Metal Inclusions on the Substrate 
Surface 


We began to present the results from our research on 
different types of metamaterials in [A47]. Then, in the paper 
[A51] we discussed the results of interesting research: two 
selected metamaterials with metal surface inclusions on a 
standard substrate have been measured by different methods 
and the extracted values of the equivalent dielectric constants 
have been compared. These materials have been developed 
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Fig. 7.1. (a) Measured relative phase delay/advance by the free-space 
method in K/Ka bands through two types of meta-samples with rectangular 
and disk metal inclusions (sample | “meta-disks”: symmetrical planar unit 
cell with disks of diameter 3 mm and distance between them 4 mm, one-side 
printed on Arlon 25N (6 mils); sample 2 “meta-surface”’: non-symmetrical 
planar unit cell 3x1.5 mm two-side printed on substrate Arlon 25N (10 
mils); (b) extracted values of the equivalent dielectric constants of the 
samples (stars — simulated by HFSS® meta-samples at 30 GHz) 
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Fig. 7.2. Split 3D models for numerical &,-extraction by free-space 
method (a) and rectangular waveguide method (b). (c) Measured TEo,; mode 
excited in a cylinder resonator with meta-samples (inset: calculated E-field 
distribution of parasitic multiple resonances around the TEo;; mode) 
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for application in antenna radomes in the Ka-band. Fig. 7.1 
presents the additional phase delay through samples’ 
thickness, measured by a free-space method in the K/Ka- 
bands, and the extracted equivalent dielectric constant £4 of 
samples, which photos have been given as insets. The first 
dielectric (“meta-disks’’) has an increased dielectric constant 
in parallel directions (measured by the two-resonator method 
in [A49]). The second sample (“meta-surface”’) placed 
horizontally, i.e. the rectangular unit cell is orientated 
perpendicularly to the electric field of the incident wave, 
shows similar behaviour. However, the same sample in a 
vertical position causes opposite phase delay (phase advance; 
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Fig. 7.1a) and the extracted parallel dielectric constant has a 
negative sign — Fig. 7.1b. The extraction has been performed 
by simulation of a specially constructed split 3D model of 
the metasamples — Fig. 7.2a. The considered samples remain 
the same dielectric behaviour by applying the rectangular 
waveguide measurement method — Fig. 7.2b, which is very 
popular for metamaterials (the extracted dielectric 
parameters have been presented in [A51]). However, the 
two-resonator method cannot be used for the characterization 
of metamaterials in all cases. Usually, the resonance 
measurements by TM mode resonators are normal, while the 
measurements by TE mode resonators sometime give 
multiple resonances instead one curve for a fixed mode. Fig. 
7.2c illustrates this effect; due to resonance excitations in 
parts of the whole sample with metal inclusions a set of 
multiple resonance curves appears (as in the measurements, 
as well as in the simulations) and the right identification now 
is impossible. In such cases, we cannot determine the meta- 
sample anisotropy for negative dielectric constants; the 
reasons have been explained in detail in [A51, A59]. 


B. Metamaterials with Controllable Filling Inclusions 


In this section, we describe a successful attempt for the 
characterization of porous metamaterials using a combina- 
tion of resonance and broadband methods [A58]. The objects 
of investigations are thin (6-100 um) anodic aluminium 
oxide (AAO) layer from anodic AlO, membranes with 
formed nano-scale air-filled pores. Such samples can be used 
as template matrices for incorporation into the pores of 
different nanostructured materials such as nanodots, nano- 
wires and nanotubes with many nowadays applications [65, 
66]. At this stage of our research, we are trying only to test 
the technology and the degree and type of inclusions are not 
well controlled. We prepared AAO samples with different 
relative purity (illustrated in Fig. 7.3; the purity depends on 
the degree of different technological contaminations). These 
initial investigations aim to determine for the first time in the 
microwave range their dielectric parameters and possible 
anisotropy and how the manufacturing technology influences 
these parameters. Table 13 presents the measured by the two- 
resonator method mean values of the parallel and perpendi- 
cular dielectric parameters for 4 selected samples with 
different purity. Definitely, samples with good purity (p1,2) 
have the behaviour of pure alumina ceramics with a set of air 
pores orientated in perpendicular directions: 1.€. Sar < Eperp ~ 
9.7-9.95. Thin samples with lossy inclusions (p3) have 
opposite behaviour &ar > &erp ~ 2.8, while both dielectric 
constants increase (Ear, Eperp > 10.4-13.5) for samples, which 
include irregular metal contaminations (p4). As we can 
conclude, in all cases the measured anisotropy is very large 
and this information can be used for fine technology 
adjustment of the porous AAO membranes. 


However, such complex materials usually may have exp- 
ressed frequency dependence of the dielectric parameters and 
anisotropy and even resonance behaviour. The applied here 
two-resonator method gives enough accurate information, 
but at fixed frequencies; to obtain informative dependencies 
in bigger frequency range we have to use another wideband 
transmission-line method (IJIC). Unfortunately, on this stage 





Fig. 7.3. a) View of anodic aluminium oxide (AAO) samples with 
decreasing purity (pl, p2 — pure or relatively pure pores; p3 — pores with 
lossy magneto-dielectric contaminations; p4 — pores with metal inclusions; 
b) microscopic side view of nanopores of diameter 35 nm 


Table 13. Dielectric parameters and anisotropy of 4 typical disk AAO 
samples with different thickness (mean values for the frequency range 7-15 
GHz) (more information in [A58]). 


Samp les j fas . par /tan Ês ; par Eperpl tan ô; > perp oe y 
S0 (p1) 29.7/42 5.96/- 9.68/- -48/- 

S1 (p2) 26/30 5.00/0.011 9.95/0.0005 -66/183 
S2 (p3) 26/18 10.50/0.016 2.78/0.0175 116/-9 
S3 (p4) 26/100 13.52/0.039 10.38/0.48 26/170 


of technology, it is very difficult to print any metal layout on 
the considered AAO samples to form a measurable transmis- 
sion line or resonator. Because of that, we applied another 
strategy for measurements — to obtain the needed information 
for the dielectric parameters, when the samples cover 
standard 50-Ohms planar lines: coplanar waveguide (CPW) 
or microstrip line (MSL) printed on a standard commercial 
substrate (in our case Arlon® 25N, 0.508-mm thick). We 
already applied this broadband method (Fig. 5.4) for a direct 
determination of equivalent dielectric constants of commer- 
cial substrates [A44], for metamaterials with metal surface 
inclusions [A51, A59], for thin nano-absorbers [A61], hexa- 
ferrites [A60] and Graphene-containing materials [A59]. 


In fact, the E fields of the dominant modes in both planar 
waveguides have mixed distribution (neither pure parallel 
nor pure perpendicular). The concrete proportions parallel- 
to-perpendicular E fields in the whole transmission lines 
(substrate and air space above) on 20-mils substrates have 
been estimated [A44] as ~38:62 in MSL and ~79:21 in CPW. 
Therefore, the covered CPW will extract predominantly 
value &, close tO & ,, While the covered MSL — a mixed 
value &g between &ar and & erp. We applied an extraction 
procedure for the determination of the equivalent dielectric 
constant &, (on the base of measured additional phase delay 
in CPW or MSL) and equivalent dielectric loss tangent 
tan ô; e) (only for lossy samples on the base of measured 
additional losses). This procedure has been developed 
applying simulations of well-constructed 3D models of the 
CPW or MSL structures (see Fig. 7.4d) with and without 
sample [A58]. Nevertheless, that the AAO samples have 
been placed on the CPW/MSL conductors under equal 
pressure (Fig. 7.4c), an unavoidable air gap g, between the 
SUT and planar line conductors appears (Fig. 6.35a) due to 
different reasons: surface roughness, cleanness, curvature, 
etc. However, this gap can be accurately determined by a 
simple way: we can preliminary measure several isotropic 
materials with known dielectric parameters (e.g. PTFE; & ~ 
2.05; Polycarbonate, €. ~ 2.78; Kapton®, &q4 ~ 3.16). In the 
considered case the obtained effective air gap values are: g, 
~ 2.5 + 0.5 um for MSL and g, ~ 17.5 + 2.5 um for CPW. 


Fig. 7.5 gives the measured additional losses and phase 
delay in CPW/MSL transmission lines caused by the 
covering AAO sample. There are two ways to extract the 
dielectric parameters of AAO samples from these data: 1) to 
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Fig. 7.4. (a) Transmission-line method: cross-section view of coplanar 
waveguide (CPW) and microstrip line (MSL) covered by SUT (E — electric 
field patterns); (b) SUT covers CPW and MSL; (c) equal pressure has been 
ensured for all AAO samples; (d) measurement setup with VNA; (e) 
Stylized 3D models of CPW and MSL and E-field distribution (red colour — 
strong fields; green — medium fields; blue — weak fields) 
tune the equivalent parameters ¢,,/tano, -7 in the 3D models 
(Fig. 7.4d) until acceptable coincidence between the 
simulated and measured phase delay and losses (last for 
lossy materials only) has been reached or 2) to introduce in 
the 3D models the pairs of parameters 6&,,,/tan6, pa and 
Eperp/tanOz perp from the resonance measurements and to 
compare the simulated and measured dependencies. In this 
example, we selected the second way, taking the dielectric 
parameters from Table 13. 


Samples S1,2 show behaviour as a pure dielectric (S1) or a 
lossy dielectric (S2) — the additional losses are relatively 
small and the phase delay increases almost linearly with 
frequency increase; the practical coincidence shows that the 
dielectric parameters from the resonance measurements are 
fully correct. However, the sample S3 has behaviour as a 
conductive material. This type of behaviour has been 
expressed with a nonlinear phase delay, which can decrease 
at higher frequencies and observation of very big losses. If 
we neglect the conductivity, the calculated losses are few dB, 
considerably smaller than the measured (not shown). For that 
reason, we have to introduce additionally to the dielectric 
parameters for this sample also an equivalent conductivity 
O-,. During the simulations, we vary the o,, values and 
compare the theoretical and measured dependencies 
simultaneously for additional losses and phases. We can 
observe a strong Oz; influence for both types of curves. The 
procedure for such a conductive sample is as follows. First, 
from the dependency for additional losses in MSL structure 
(Fig. 7.55) we can evaluate the equivalent conductivity 
Oeq_msL (the dielectric constant doesn’t influence so much the 
losses). The obtained values are: O4 ms ~3-10 S/m (for the 
frequency range 1-5 GHz), ~10-20 S/m (5-10 GHz), ~20-35 
S/m (10-30 GHz), ~35-50 S/m (30-40 GHz). The simulated 
values for Oz, cepw in CPW structure are similar (Fig. 7.5a); 
typically O4 cew ~35 S/m. Now, if we consider the phase 
delay dependence in MSL, the combination of values &¢ ms 
~9-10 and oy, usz ~35-40 S/m is the solution for this con- 
ductive sample S3. The phase delay dependence in CPW 
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Fig. 7.5. Measured additional losses and phase delay in CPW (a) and 
MSL (b) transmission lines after covering by thicker AOO samples (solid 
curves). Theoretical dependencies (dashed curves), obtained by 3D 
simulations with dielectric parameters from Table 13 (anisotropy option) 
and setting of effective air gaps g, = 2.5 um (MSL) and 17.5 um (CPW) 


gives slightly different combination, &, cpw ~11-12 and 
Oeg_cpw ~50 S/m; the difference means that the parallel 
dielectric constant and conductivity are bigger. 


The obtained results for the dielectric characterization of 
AAO samples show that the proposed combination of 
resonance and broadband microwave methods is working 
well for such complex metasamples with variety of dielectric 
or metal inclusions. Samples with a small amount of metallic 
contamination can be better characterized with high 
accuracy. Samples with higher conductivity are also 
measurable: dielectric constants and anisotropy by the 
resonance method and equivalent conductivity — by the 
broad-band transmission-line method. All these interesting 
and promising investigations are in progress. 
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C. Carbon-Containing Materials 


A good example for successful simultaneous determinati- 
on of the pair of equivalent parameters &,/o., according to 
the Drude’s model for conductive dielectrics (when the loss 
tangent has been neglected; expression (6)) 1s the character1- 
zation of several carbon-containing materials. In this section, 
we Selected three types of samples: spectral pure Graphite, 
air-filled Graphite and pressed Graphene. The measured 
additional phase delay in CPW transmission line has been 
given in Fig. 7.6a, while the extracted equivalent dielectric 
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Fig. 7.6. (a) Measured relative phase delay/advance in Graphite, air-filled 
Graphite and pressed Graphene samples by covered CPW (sample diameter 
~9 mm; thickness 0.9 mm; g, = 20 um gap “sample-conductor”’; (b) 
extracted &, values (inset: dielectric constant of graphene by [66]) 
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Fig. 7.7. (a) Measured additional losses in CPW with carbon-containing 
samples as an overlay; (b) extracted frequency dependencies of &,; stars: 
values of the equivalent dielectric constants by taking into account of the 
actual conductivity (inset: extracted equivalent conductivity o,) 


constant &, — in Fig. 7.6b (when O has been neglected) 
[A51, A59]. The air-filled Graphite has the behaviour like a 
low-Dk dielectric; é,~1.4-1.5). However, the other materials 
demonstrate completely different behaviour: the phase delay 
at lower frequencies transforms into phase advance (positive 
delay). In this case, the extracted ¢,, values show strong fre- 
quency dependence — from 15-20 below 2 GHz up to | and 
even values below | for both Graphite and pressed Graphene 
samples. This dependence corresponds to the published 
results from other authors (e.g. [66]). However, the measured 
additional losses in CPW of these samples are very big (Fig. 
7.7a) and cannot be explained with dielectric losses only; the 
samples have also a big equivalent conductivity o.,. As in 
the previous example, we can extract this parameter (the 
differences between samples are well observable; we added 
also the dependence for Nitrogen (N)-doped Graphene [67]; 
5.5-% doping level). We established that the sensitivity and 
accuracy of the CPW measurement systems to the Ox 
variations are better in comparison to the MSL case; that's 
why we concentrate our efforts on CPW structure only. The 
extraction procedure is as follows. We use in the 3D models 
(inset in Fig. 7.7b) the initial ¢,, value (determined for oj, = 
0; Fig. 7.6b) and starting to tune the conductivity to match 
the simulated and measured losses at a fixed frequency, we 
can determine o,, at the same frequency. Then, simulating 
the phase delay/advance with the new © value at each 
selected frequency, we can recalculate the final ¢, value — 
results for the extracted &,/o., dependencies are visible from 
Fig. 7.7b (the changes for &4 are negligible). The measured 
Oe, Values in the frequency range 5-40 GHz are relatively 
small due to the air gaps in the weakly-pressed samples. 


D. Plasmas 


Plasmas (gaseous, solid-state) are such media, which also 
can have an expressed anisotropy of their dielectric 
properties depending on the realized type of plasma 
discharge and orientation of the external magnetic fields (if 
exists). In our research, we developed several hairpin probes 
for resonance measurements of plasma densities in different 
directions [A27, A37, A37]. The hairpin probe is a simple 
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Fig. 7.8. (a) U-shape straight (180-deg) hairpin probe and (b) its 3D 
model in an eigenmode option of the simulator; (c) U-shaped hairpin probe 
orientated at 90 deg with ceramic support on the shorted end; (d) measured 
return losses of the hairpin probe in the air with the exciting first three 
TEM-mode resonances (2-10 GHz) 
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Fig. 7.9. a) Tools for mounting of the hairpin probe in the vacuumed 


plasma reactor; c) dependence of the hairpin resonance frequency fhaipin ON 
the plasma density ne. Inset: expression between fhaipin and ne 
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Fig. 7.10. (a) Dependence of the hairpin resonance frequency fhaipin on the 
effective plasma permittivity & and simulated E-field distributions in 
different points. Inset: expression between f, and ne; (b) the same for the first 
two modes in the hairpin probe (Inset: expressions for the relations “sp — 


Shaipins fr" 

quarter wavelength two-wire resonator, which can be relati- 
vely easy designed for TEM-mode operation [68, A27]. The 
open end of the structure ensures a maximum of the E field 
of the standing wave, which makes this resonance probe 
enough sensitive to the changes of the dielectric parameters 
of surrounding medium (e.g. effective plasma permittivity 
Ep ef). Lhe opposite short end has a maximum of the H field, 
which allows to achieve a quite stable and reliable magnetic 
coupling between the resonator and the feeding coaxial cable 
using H-type coaxial loop probe — see Fig. 7.8. The main 
benefit of this method in comparison with the other methods 
for plasma density evaluation is the fact that the measure- 
ments are based on a determination of the resonance 
frequency fhaipin instead of any absolute measurements of E- 
filed magnitudes, which ensure better sensitivity, accuracy 
and measurement simplicity. 
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Fig. 7.11. Electron plasma density n, versus the reactor power (W) of 
100-% oxygen (Oz) plasma (a) and of 100-% SF, plasma (b) at constant 
pressure (100 mTorr) 

Fig. 7.9c illustrates the dependence of the hairpin 
resonance frequency fhaipin on the electron plasma density ne 
in a first-order approximation; similar dependence allows 
after calibration on air the relatively accurate determination 
of ne during different plasma technological processes in 
vacuumed plasma reactors. The reason for existing of this 
dependence is the relation between the effective plasma 
permittivity ¢ , and of the ratio between the plasma 
frequency f, and the hairpin resonance frequency fhaipin (See 
the inset in Fig. 7.10). In our research, we managed to 
determine the zone for reliable usage of this dependence (in 
the interval &p c (1, —1); see Fig. 7.10a) between the area of 
saturation and the fact that the sensitivity of the hairpin probe 
increases by using of the high-order modes (Fig. 7.10b), 1f 
they have been well identified [A35]. 


We developed an experimental hairpin probe [A27] for 
implementation in plasma reactors with aggressive gases (as 
in [69]) like oxygen O, and SF. The developed probes have 
been incorporated in some reactors of Oxford Instruments© 
(see [70]). The special features of this type of hairpin probes 
are the type of used probe metal (Platinum + 30% Rhodium), 
which ensures bigger sensitivities than the conventional 
metals (see dependencies in Fig. 7.11) and the feeding 
coaxial cable filled with S102, which is working at relatively 
high temperatures in the considered plasma reactors. 


E. Experimental Observation of Anisotropic Magneto- 
Dielectric Effects in a Class of Artificial Materials 


In IIC, we introduced the material properties of the so- 
called magneto-electric materials and the possibility for their 
characterization. These specific materials will pay increasing 
role as integrated tunable components in the modern 
microwave devices [71] due to the ability to separately 
control their magnetic properties by an external dc voltage 
and their electrical properties by an external magnetic field. 
However, the experimental separation of these effects is not 
an easy task. In our research, we already succeeded to detect 
magneto-dielectric effects in several artificial materials in the 
frequency range 2-40 GHz [A59, A63, A68]. In this paper, 
we will briefly describe as an example the experimental 
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Fig. 7.12. Applied external dc magnetization of sample in parallel/ 
horizontal (a) and in perpendicular/vertical (b) direction; (c) electromagnet 
system. Legend: 1 — sample; 2 CPW substrate; 3- electromagnet 
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Fig. 7.13. Measured dependencies for the additional phase in CPW 
structure with multiferoic sample in vertical (a) and horizontal (b) magnetic 
field Hp = 0-1.33 kOe, normalized to values at Ho = O (the red arrows show 
the directions for Ho increasing). 
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Fig. 7.14. Resonance curves of different modes in rectangular resonators 
(WR28) with prisms (1.22x2.07x10.16) from multi-feroic in the centre, 
which have been influenced only by the dielectric (eps) (a) or magnetic (mu) 
(b) constants of a sample in horizontal magnetic field Ho 0-1.33 kOe (the red 
arrows show the directions for Ho increasing). 
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Fig. 7.15. a) Variations of the additional phase in CPW structure with 
multi-feroic sample at 30 GHz v/s magnetic field Ho; b) shift of resonance 
frequency of 6 TEjo, modes due to magnetized multi-feroic prism in external 
magnetic field Hp = 0-1.33 kOe, normalized to values at Ho = 0 
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verification of supposed magneto-dielectric effect in a 
multiferroic material Sr3Co.Fe2404, ([A60], see also VB). 


In most cases, the magneto-dielectric effect in multiferroics 
has been measured at relatively low frequencies, (typically 
several hundred MHz, [72]). Only a few papers consider the 
characterization of this effect in the microwave range [73]. 
Here we describe an experimental characterization of 
dielectric/magnetic constants’ variations of Z-type hexa- 
ferrite Sr3Co.Fe24O,4; in an external magnetic field. Fig. 7.12 
shows the applied setup that combines measurements in the 
frequency range 0.01-40 GHz of CPW transmission line 
covered by magnetized sample by parallel or perpendicular 
external magnetic field Hj between O up to 2.4 kOe. The 


samples have been directly placed on the CPW conductor 
and the additional losses and phases caused by the 
magnetized sample can be measured, when A increases. Fig. 
7.13 a,b illustrates the differences in additional phase for Ho 
in both horizontal/vertical directions. The influence of the 
horizontal Hp is 5-6 times bigger, which clearly indicates 
anisotropy of the measured effect — Fig. 7.15a. However, we 
cannot separate in these measurements the dc magnetic field 
influence on the sample dielectric and magnetic constants. 
This information we can obtain by realized resonance 
perturbation measurements in a rectangular resonator with 
small ferrite prisms in the resonator centre. From the changes 
in the resonance frequencies of the excited TE), modes (Fig. 
7.14), we can separately extract the variation of the equiva- 
lent magnetic Leg (P = 2, 4, 6) and dielectric &, (p = 1, 3, 5) 
constants of the ferroic sample. The analysis of the results in 
Fig. 7.15b allows us to conclude that the smaller external 
magnetic fields (Hp < 0.4 kOe) influence stronger the dielec- 
tric constant &€,, while the stronger fields (Hp > 0.8 kOe) — 
the magnetic constant Leg, Where £g variations show a satu- 
ration. These conclusions have been confirmed by the broad- 
band CPW measurements; see dependencies in Fig. 7.15a, 
where a 1.4-8.5 % variation of the material constants can be 
achieved in the frequency range 28-38 GHz, which is enough 
for possible applications as controllable delay lines. The de- 
tailed investigations of the magneto-dielectric effects in fer- 
roic and other types of materials are still in progress [A68]. 


VU. EQUIVALENT CONDUCTIVITY OF MATERIALS 


Nowadays the conductivity of materials starts to play an 
increasing role in the design of new RF devices, comparable 
with the importance of the role of the actual dielectric and 
magnetic parameters in the same design. The microwave 
engineers need to know the actual equivalent conductivity of 
the metals and metallization layers used in different bulk and 
planar devices accurately to determine the conductor losses 
and reflections. Moreover, today many devices are realizing 
with metalized injection-moulded or 3D printed dielectrics 
instead solid metal walls, where the equivalent conductivity 
depends not only on the used metal but also on the surface 
roughness, cleanness, used protective layers, etc. Many new 
materials in modern electronics have smaller than ordinary 
metals but noticeable conductivity — carbon-containing 
materials, metasurfaces, absorbers, conductive textile fabrics, 
etc. In this section, we present our experience in the area of 
characterization of equivalent conductivity in the interval 
from typical values for conductive metals (e.g. ~4-6.10’ S/m) 
up to 10°-10* S/m by resonance methods and up to 1-10 S/m 
by broadband methods. 


A. Equivalent Conductivity of Metallized Plastic Samples 


The nowadays technologies for creating of innovative 
antenna components for airborne, satellite and different 5G 
applications — metallized plastic details, 3D printed 
structures, metamaterials with metal inserts, metallized 
composites or textiles, etc. sharply increased the needs to 
characterize the actual conductivity of these materials due to 
the aims for more reliable 3D design. Actually, the deter- 
mination of the conductivity Oyw in the microwave range 
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Fig. 8.1. (a) Resonators with decreasing “volume-to-bottom surface” 
ratio; excited TE;o5 mode in rectangular Rs1 (b) and sphere-cylinder Rs2 (c) 
measurement resonators; (d) reference curves “o.,-Qo* with Brass reference 
for determination of equivalent conductivity of surface under test (SuUT) 
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Fig. 8.2. Calibration curves “equivalent conductivity — unloaded Q 
factor” for the all-Brass resonator and “Brass conductivity reference” (a) and 
for the all-Silver resonator and “Silver conductivity reference” (b) and test 
values of conductivity for several known conductive materials. 

(instead of the widespread dc Oge or optical oœ values) is a 
classical problem. It is also a known fact that the so-called 
effective conductivity diminishes for electroplated metal 
surfaces [74]. It depends not only on the skin depth but 
additionally on the surface cleanness, roughness, flatness, 
coatings, technological scratches and other surface 
irregularities. For this reason, we developed in [A43] an 
efficient microwave resonance method for determination of 
the useful parameter equivalent conductivity Oeg, which 
better characterizes the metal surfaces in X-Ka bands and 
more or less differs from the dc bulk conductivity usually 
introduced in the 3D simulators for antenna components’ 


design. The idea of the method is based on the replacement 
of one flat surface of a selected volume resonator with the 
surface under test (SuUT). To increase the resonator 


sensitivity to the conductivity variation of SuUT, we started 


to search specific measuring resonators with small “volume- 


to-bottom 
resonators — the standard rectangular resonator with excited 


surface” ratio (Fig. 8.la). We selected two 
TEiop modes (p = 1-6) (Fig. 8.1b) and a sphere-cylindrical 
resonator with excited similar modes (Fig. 8.1c). Applying a 
“Brass conductivity reference” for calibration of the 
measurement systems, we managed to determine with 
satisfying accuracy the actual equivalent conductivity Oeg Of 
many metalized surfaces (e.g. Alodine 1200 passivation 
layers over Dural surfaces; Fig. 8.1d) used in different 
antennas. However, the benefit of using the sphere-cylinder 
measurement resonator could be minimized due to the 
complex E field distribution for higher modes and the 
“detachment” of that field from the bottom resonator wall 
(SuUT). Actually, the applying of a “Silver conductivity re- 
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Fig. 8.3. (a) Measuring resonator with SuUT; (b) Two measuring 
resonators Rs1 and Rs2; (c) Reference surfaces from Al, Cu and Ag; (d) 
Surface view of several samples: Reference Silver EP; bulk Ag; JM — 
sprayed Ag nano-coating; 3D EP — electro-plated Ag on 3D printed 
photopolymer; EL2 — electro-less Ag; EL2+nano protective layer; EL2+Rd 
Rhodium protective layer; (e) smooth and rough Al surface 









































7.0 
can 6.5 
= Lal 
= 604° 
“A . 
~a J. 
= 5.0 
b 4.5 
4.0 
aio 
3.0 
ss | Lit 
14 16 18 20 22 24 26 28 30 32 34 36 38 40 (a) 
f . GHz 
7.0 7.0 
6.5 6.5 = 
— 60F 46.0 5 
& 5s dss 8 
Yi 5.0 5.0 y 
o 45E Je J45 g 
a r 3.5 3 
b” 3.0 — @ = m 430 6 
2.5 Attenuation in WR28 2.5 Q 
2.0 Equivalent conductivity a at 32 GHz i 2.0 Zz 
1.5 - — able level = | 0dB m z — $ — 1.5 
1.0 < 1.0 
0.5 F 0.5 





m, 
_— 
=- 
—_ 


0.0 (b) 
samples 
Fig. 8.4. (a) Measured equivalent conductivity 6., of Ag-coated samples 
by different technologies (see the legend in Fig. 8.3d; all Ag layers are 5-12 
um thick); (b) Mean o,, (at 32 GHz) and calculated attenuation in WR28 
waveguide for Ka-band for the considered samples 
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reference” for the measurement resonators can increase the Q 
factors of the excited modes twice and to additionally 
increase the sensitivity [A55]. The older calibration curve by 
the “Brass conductivity reference” covers the interval for 
accurate 6,, determination from 7.10’ to 5.10° S/m, while the 
new calibration curve by “Silver conductivity reference” 
covers considerable bigger interval — from 7.10’ up to 10° 
S/m — see the set of calibration curves in Fig. 8.2. The 
procedure for the determination of equivalent conductivity 
6-, of SUUT is as follows. First, the unloaded Q factors of all 
excited TE modes have to be measured in the resonator with 
silver-plated walls, including the “Silver reference” as a 
bottom wall (instead of SuUT), Fig.8.3a,b,c. The simulations 
of the 3D models for each mode (Fig. 8.1b,c) allow 
determining the equivalent conductivity of all resonator 
walls. Then the “Silver reference” has to be replaced with the 
corresponding SuUT and the new unloaded Q factors to be 
measured. The new simulations by varying the conductivity 
only for the SuUT allow to achieve practical coincidence 
between the measured and simulated Q factors and to obtain 
6-, Of SUUT in the considered frequency range. 


Applying this method, we succeeded to characterize the 
actual conductivity Oe of a lot of artificial metalized 
surfaces. First of all, we investigated the influence of the 
surface roughness, described by the parameter R, — the total 
“peak-to-valley height” of the surface bulges. Fig. 8.3e 
shows two Dural surfaces: 1) “smooth” (with R, < 0.05 um) 
and 2) “rough” (with R, ~0.8-1.6 um). We established a 
decrease of o,, with ~10 % in Ku-band and ~25 % in Ka- 
band of rough surfaces compared with the smooth surfaces 
of the same metal sample [A41], which proves the acceptable 
sensitivity of the proposed measurement method. 


Fig. 8.3d illustrates the view of the top surfaces of several 
silver-coated plastic samples, prepared by different 
technology; o,, dependencies have been drawn in Fig. 8.4a. 
Acceptable values Oe = 4.6-6.2 S/m are fully realizable by 
electro-less (EL), electro-plated (EP) or nano-Ag-aerosols- 
sprayed (JM [76]) technologies. The fine surface finishing 
(EL1) increases the conductivity with 15-20% (EL2); the 
protective layers decrease o,, with ~10 % (EL1 + Rhodium 
protective layer) or up to 35 % (ELI + protective nano- 
coating). However, the importance of protective layers 
becomes visible in the environmental tests; an 85 % decrease 
of the equivalent conductivity in Ku-band was detected for 
pure electroless silver surface (EL1) in a salt-frog test (from 
6.2 to 1.2 S/m [A55]), while the protected with Rd surface 
keeps value o,, ~4.3 S/m. Similar measurements allow us to 
establish the acceptable equivalent conductivity for 
metalized injection-moulded or 3D printed plastics with 
metallization; the calculated attenuation ~1 dB/m in WR28 
waveguide for Ka-band can be satisfied for mean Ge 22.5 
S/m (at 32 GHz) — see the dependencies in Fig. 8.4b. 


B. Equivalent Conductivity of Both Sides of Metallization 
used for Reinforced Commercial Substrates 


The determination of equivalent conductivity of copper 
metallization used for the microwave substrates is extremely 
important especially because PCB producers apply Cu folio 
with different roughness of both sides [76]. In the considered 


in the previous section waveguide components made from 
metalized plastics, the losses depend on the equivalent 
conductivity Oe of the metallization top surface. 
Contrariwise, this role for the substrate metallization 1s 
playing from the bottom side (to the substrate) and strongly 
determines the transmission-line losses of printed planar 
structures, especially in the mm-wavelength range [77]. To 
minimize these losses, but to ensure acceptable adhesion of 
the metallization to the substrate surface, the producers apply 
different techniques to increase the roughness of the bottom 
surface, which decrease the equivalent conductivity. Our 
method [A55] allows accurate distinguishing of this 
difference — see the measured dependencies of the top and 
bottom sides’ conductivity o,, for several realistic Cu folios 
in Fig. 8.5a; the top surface has typically 3-5 times higher 
Ovg, excepting reverse treated substrates. The reasons are well 
explainable. The standard technology for electro-deposited 
(ED) Cu folio ensures typically RMS irregularity Rrys ~0.3- 
0.4 um (which better represents the surface roughness RGH) 
of the metal top surface (drum or resists side), while the 
values for the opposite backside are Rrys ~ 1-4 um [76]. 
Only by R/A (rolled/annealed) techniques, these values can 
be decreased to Rrys ~0.4 um. The producers verified the 
obtained volume and surface resistivity by the IPC-TM-650 
2.5.17.1 test method [78] at low frequencies. Our 
measurement method allows more reliable characterization 
of the equivalent conductivity of both metallization sides at 
the actual working frequencies. Fig. 8.5b shows the 
measured attenuation of two commercial substrates with Dk 
~3 with ED and LoPro® metallization and calculated values 
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Fig. 8.5. a) Measured equivalent conductivity of copper metallization for 
microwave substrates; top/resist (solid curves) and bottom (dashed curves) 
sides (17.5-um thick Cu folio). Legend: 1 — Rolled/Annealed (R/A, Rrums = 
0.4 um); 2 — Reverse Treated (RTF); 3, 4 — Electro-Deposited ED with 
different adhesive glue; 5 — ED, additionally electro-plated (35-um thick). 
The pictures present the images of the bottom surfaces (to the substrate); b) 
Attenuation in 50-Ohms microstrip line on substrates with ED or LoPro® 
metallization, calculation by measured o,, (see Table 14) 
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Table 14. Measured equivalent conductivity o., of the Copper folio for 
microwave substrates and calculated conductor and total losses in 50-Ohms 
MSL (insets: RGH of both metallization sides). 





Meisl type RGH Rao) Sim aces ai dBiem 
Ideal Cu 5.80.107 0.041 0.065 
R/A (0.4 um) 5.09.10 0.044 0.067 = Rays ~ 0.3- HII 
ED LoPro (0.9m)  1.37.107 0.084 0.112 0.4 um 
ED (2.0 um) 0.63.107 0.113 0.161 


Substrate side > MASNA 


RTF (1.8 um) 0.79.10° 0.110 0.138 





Suid Sli 
at 32 GHz for the same cases by TRL calculators. The 


rolled/annealed (R/A) metallization used for substrates in the 
mm-wavelength range has better o,, (~1.5x10’ S/m) in 
comparison to the electro-deposited metallization (~0.7x10' 


S/m) (Fig. 8.5a), which ensure attenuation decrease from 
0.16 dB/cm to ~0.11 dB/cm (Fig. 8.55; Table 14). 


C. Equivalent Conductivity of Meta-Surfaces and Materials 
with Low Conductivity 


The proposed in [A55] method for accurate 
characterization of equivalent conductivity even for values 
less than o., < 10°+10' S/m is fully applicable for different 
metasurfaces and materials with low conductivity. A good 
test of the used measurement resonators is the 
characterization of carbon-based samples. Fig. 8.6a shows 
how the measurement resonators with high Q factors allow 
reliable o,,-measurements of carbon-based samples (ordinary 
or carbon nanotube CNT cloths with epoxy; the CNT 
samples have 10 times higher Oe ~ (0.2-1)x10° S/m instead 
of ~(0.4-1)x10* S/m. The conductivity of the textile fabrics 
with different conductive fibres content is also well 
determinable —10°-5.10° S/m. Very successful could be the 
Og Characterization of meta-surfaces — a possible resonance/ 
periodical conductivity behaviour is fully detectable and 
measurable for the already considered metasurface samples 
(compare with the dependencies for the resonance behaviour 
of dielectric constant in Fig. 7.15). 
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Fig. 8.6. a) Measured equivalent conductivity of low-conductive samples. 
Legend: 1 — pure Carbon cloth; 2 — Carbon cloth + epoxy (1:1); 3 — CNT 


cloth + epoxy; 4, 5 — conductive textile fabrics (90% and 50 % conductive 
fiber content); b) 6 — artificial metasurface; 7 — meta disks 





IX. CONCLUSIONS 


In this paper, we summarize the main results from our 
many-years’ research work devoted to the characterization of 
the dielectric properties of materials used in modern 
microwave electronics. In the first part, we prove the 
assumption that the well-determined complex dielectric and 
magnetic constants of the materials are important not only 
for improving the accuracy of the modern 3D design of 
nowadays electro-dynamic structures but also because the 
integral character these parameters ensure valuable 
information for the real material composition, structure, 
character of inclusions, building blocks and unit cells 
orientation, used technology, conditions for the material 
preparation, etc. We have shown that additional very useful 
information can be achieved when the actual anisotropy of 
the material constants has been determined and compared — 
different behaviour of their permittivity/permeability in 
different directions. In many _ standard applications 
(excepting usage of some metamaterials), the anisotropy of 
the most of artificial materials is more or less undesired 
properties, which decreases the accuracy for the electro- 
dynamic description of these materials (e.g. in the cases of 
reinforced substrates, wearable textile fabrics, antenna 
radomes, 3D printed dielectrics, etc.). However, it is 
important for us that nowadays most of the PCB 
manufacturers, EM simulators’ developers and most of the 
material users realize the importance of the actual anisotropy 
of the modern materials and all of them try to solve and to 
take into account this property. 


Our numerical models and experimental methods for 
characterization of the sample dielectric and magnetic 
properties including their anisotropy give satisfactory 
accuracy practically for all possible applications. The 
relatively full set of the implemented resonance and broad- 
band measurement methods in our Microwave laboratory 
gives us the possibility for deep investigation and 
characterization of a variety of different materials (some of 
them completely new) in the nowadays electronics — 
microwave reinforced substrates, ceramics, multilayer 
composites, different dielectric mixtures, 3D printed 
dielectrics, textile fabrics, metamaterials, thin micro- and 
nano-films, carbon-containing materials, fresh and dry plant 
tissues, plasmas, etc. 
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